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^ Abstract 

•i-H 

j_j ■ The core of this article is a general theorem with a large number of 

^ specializations. Given a manifold and a finite number of one-parameter 

groups of point transformations on A^ with generators F, ■ ■ ■ we 
obtain, via functional integration over spaces of pointed paths on A^ (paths 
with one fixed point), a one-parameter group of functional operators acting 
on tensor or spinor fields on A^. The generator of this group is a quadratic 
form in the Lie derivatives C.X(a) X(Q,)-direction plus a term linear in 

The basic functional integral is over L^'^ paths x : T ^ A^ (continuous 
paths with square integrable first derivative). Although the integrator is in- 
variant under time translation, the integral is powerful enough to be used for 
systems which are not time translation invariant. We give seven non trivial 
applications of the basic formula, and we compute its semiclassical expansion. 

The methods of proof are rigorous and combine Albeverio H0egh-Krohn 
oscillatory integrals with Elworthy's parametrization of paths in a curved 
space. Unlike other approaches we solve Schrodinger type equations directly, 
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rather than solving first diffusion equations and then using analytic continu- 
ation. 
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I - Introduction 



We have studied many applications of functional integration looking for 
its substantifique moelle^ . Little by little, several ideas have taken shape 
concerning the domains of integration, the integrators, and the integrands. 

1. The domain of integration is a function space. 

Working with an infinite-dimensional space is easier than working with 
the limit for large n of the product of n copies of a finite-dimensional space. 
For example, a space of pointed paths (paths with one fixed point) is con- 
tractible even when the paths take their values in a non contractible space. 
Over the years the advantages - often the necessity - of working with spaces 
of paths rather than with the discretized version of the paths have become 
increasingly apparent: 

- Semiclassical approximations, even in the presence of caustics, are ob- 
tained by expanding a functional on the space of paths around a dominating 
contribution [1,2]. 

- If the paths take their values in a multiply-connected space the topology 
of the space of paths plays a central role [3] . 

- A change of variable of integration, regarded as a map on the domain of 
integration gives, in two lines, the Cameron-Martin formula [4] obtained from 
discretized paths via a lengthy derivation, and generalizes its applications [5] . 

- Computing functional determinants using properties of linear maps on 
Banach spaces is simpler than computing limits of finite-dimensional deter- 
minants [6]. 

- A major progress in the definition and computation of functional inte- 
grals was achieved by Elworthy [7] when he parametrized the space of paths 
on a Riemannian manifold M, with fixed initial point a, by the space of paths 
on TaM starting at the origin of T^M. 

The importance of the domain of integration, noted in the above exam- 
ples, is even more striking in the formulation presented here. A key point is a 
generalization of Elworthy 's idea. Consider a finite-dimensional manifold N, 
and denote by T some finite time interval. Here we parametrize a space Vx^N 
of pointed paths a; : T — > A'' by a space VqW^ of pointed paths ^ : T ^ M^. A 
general construction for maps 

Vo^^ ^ 7^0=0^^ by z^x 

^ F. Rabelais. Literally "bone marrow as a producer of substance" in 
Gargantua, Prologue de I'auteur. 
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is given in Section II by solving suitable first-order differential equations, not 
by constructing stochastic processes. It is stated in terms of vector fields Y, 
• • • , X(^(i^ on A^. By specializing N and the vector fields, one basic func- 
tional integral yields, with rigor and no Ansatz, a great variety of functional 
integrals which are solutions of complex problems. A number of them are 
presented in Section IV. 

In general, to define the domain of integration, one needs to specify: 

- the analytic nature of the paths^ (continuous, L^, L^'-^, • • •); 

- the domains of the paths and their ranges; 

- the behaviour of the paths at the boundaries of their domains. 

2. Integrators cannot be expected to be universal. 

The naive approach to the definition of an integral in an infinite number 
of variables is to take a limit d = oo in a d-dimensional integral. Because of 
scaling problems this procedure is well known to abort. For instance, if we 
wish to evaluate (for a > 0) the integral 

loo •■= [ d°°xexp(--\x\'') , (I.l) 
Jroo \ a / 

oo 

where |xp := (^")^; we may first evaluate the corresponding integral 

a=l 

Id := / dx e^p ( --\x\'^) (1.2) 
jRd V a / 

for finite d and then set d = oo. Since 1^ = a^^"^ we get 

C if < a < 1 
/oo = < 1 if a = 1 (1.3) 
[ oo if 1 < a 

and this fails to be continuous in the parameter a, as should be reasonably 
desired. 

A way out of this diflficulty is to introduce, for each value of the scaling 
parameter a > 0, an integrator VaX in the d-dimensional space W^, namely 

VaX = a-'^/^dx^---dx'^ (1.4) 

2 A path ^ : T ^ is said to be L^'i if dt |i(t)|2 < oo where z{t) = 
dz{t)/dt. 
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and to remark that it is characterized by the following integration formula 
/ T>aX ■ exp (— — |a;p — 27rz {x',x)j = exp (— Trajx'p) . (1.5) 



Here x' runs over the space dual to M*^ and the scalar product is given^ by 

\ t/^ ^ tXy j tt/ tXj • 

This formula is dimension-independent and hence suitable for the gener- 
alization from M'^ to a (real) Banach space X. Let X' be its dual and consider 
two continuous quadratic forms Q on X and on X'. Assume that Q and 
W are inverse to each other in the following sense. There exist continuous 
linear maps 

i:> : X ^ X' , G : X' ^ X 

such thaf^ 

{DG = GD = 1 
{Dx,y) = {Dy,x) (1.6) 
Q{x) = {Dx,x) , W{x') = {x',Gx') . 

Here {x', x) denotes the duality between X (elements x) and its dual X' 
(elements x'). 

Then we define the integrator Vg^QX (also denoted Vx for simplicity) by 
the following requirement 

J Vs^QX-exp(^--Q{x)-27ri{x',x)^ =exp{-TvsW{x')) . (1.7) 



Here x' runs over X' and there are two cases: 

i) s = 1 and Q is positive definite, namely Q{x) > for a; 0; 

ii) s = z and there is no restriction on Q except it be real. 

Introduction of the parameter s enables us to treat in a unified way the dif- 
fusion and the Schrddinger equations^. We can introduce a suitable space 

^ Here and in the rest of this paper we use the Einstein summation con- 
vention over repeated indices. 

^ In standard applications, Z) is a differential operator and G the corre- 
sponding Green operator, taking into account the boundary conditions of the 
domain of D. 

^ It is not true, as is still often stated, that the case s — i has no mathe- 
matical foundation. See for example references [8, 9, 10, 11]. 
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.F(X) of functionals on X integrable by "Ps.Q and a norm on .F(X), and then 
compute integrals of the type 

/ = / V.^QX ■ F{x) (1.8) 

for F in J^(X). In both cases, s — 1 and s = i, we shall call Vg^QX a Gaussian 
integrator. 

Gaussian integrators have the following properties: 

T>{x + xq) = T>x, xq a fixed element of X (1.9) 

V{Lx) = I Det L| • Pa;, L : X ^ X (I.IO) 

where L is in a suitable class of linear changes of variables of integration, 
including the obvious case where Q{Lx) = Q{x) and | DetL| = 1 (cf. formula 
(A. 140) in Appendix A). 

Gaussian integrators are not the only possible integrators. In reference 
[11] we have developed an axiomatic for functional integrals on a Banach space 
$ expressed in terms of integrators Vq^z defined by 

/ Pe,zv^-e(^,J) = Z(J) (1.11) 

for (/9 in $, J in the dual $' of $, where Q and Z are two given continuous 
bounded functionals 

: $ X ^ C , Z : ^ C. 

In quantum field theory, we interpret (/9 as a field and J as a source, Z{J) is 
then the Schwinger generating functional for the n-point functions. 

3. Integrands and integrators. 

Splitting the quantity inside the integral sign into "integrator" and "in- 
tegrand" belongs to the art of integration, but rules of thumb apply: 

- When the functional integral has its origin in physics try not to break 
up the action into, say, kinetic and potential contributions. On the other 
hand, do not hesitate to work with a potential which is a functional of a path 
rather than a function of its value (e.g. in equation (III.l), V is a functional 
of z, not a function of z{t)). 

- Look for a possible change of variable of integration; this may suggest 
a practical choice for the integrand. 
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- Gaussian integrators have a wealth of simple, powerful properties; look 
for exponentials of quadratic forms, this suggests a practical choice for the 
integrator. 

4. A basic functional integral. 

The core of this article is a theorem which provides the mathematical 
underpinning for a great variety of functional integrals. It consists of two 
parts: the definition of a functional integral, and the partial differential equa- 
tion satisfied by the value of the functional integral, as a function of a set of 
parameters. Given a manifold consider the L^'^ paths over a finite time 
interval T with values in N, x : T ^ N, and with a fixed point xq G N; for 
instance, if T = [taj^b]; the fixed point can be either x{ta) or x{tb). Given 
d + 1 vector fields Y and X(^a) on N, define a map P from a space ■Po^^'' of 
L^'-*^ paths into a space Vxo^ of L^'-*^ paths, 

P : VoW" ^ Vx,N, 

hy P{z) = X, where 

dx{t, z) = z))dz'^ + Y{x{t, z))dt 

x{to,z) = Xq. 

Here to and xq are fixed, with to in T and xq in A^, and the paths ;2 : T — M'^ 
satisfy z{to) = 0. In general, the vector fields do not commute, that is: 

[^(a),^(/3)] 7^0 , [r,X(,)]^0; 

therefore the solution of (1.12) is of the form 

x{t, z) = XQ ■ J:{t, z) (1.13) 

where x is a function of xq and t, and a functional of z. Here T^it^z) 
is a transformation in N ^ depending on t and z as stated. Only when 
= [y, = can one express S(t, z) as a function of t and 

z{t). 

In Appendix C the familiar theorems of differential equations are ex- 
tended to cover the properties of (1.12). 

Given (1.12) one can express an integral over Vxo^ ^ integral over 
VoM.'^, which is an integral that one can manipulate and compute. The partial 
differential equation satisfied by the integral on Vxt^N is expressed in terms 
of Lie derivatives along the vector fields Y and Conversely, given a 

parabolic partial differential equation, one can construct in many cases the 
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path integral representation of its solutions (see for instance in Section IV, 
the lift of a covariant Laplacian at a point of a Riemannian manifold in terms 
of Lie derivatives at a point of its frame bundle [13]). 

The basic equations are given in the first paragraph of Section II, followed 
by a summary of notations used throughout the paper. 

5. Examples. 

In Section III we compute semiclassical approximations of the basic func- 
tional integral (11. 1), and in Section IV we specialize the basic integral by 
making particular choices of the manifold N. We treat in detail the following 
examples: 

- TV = M'^, in cartesian or polar coordinates. 

- is a frame bundle 0{M) over a Riemannian manifold M. We give 
the explicit functional integral representing the solution of the Schrodinger 
equation on M, with initial wave function 0. 

- N is a multiply-connected space. 

- N is a t/(l)-bundle; the basic integral solves the Schrodinger equation 
for a particle in an electromagnetic field. 

- N is a symplectic manifold; coherent-state transitions can be obtained 
from the basic integral. 

We conclude this section by an analysis of the Bohm-Aharonov effect, where 
all the previous techniques are brought to bear. 

6. Techniques. 

In the course of computing our basic functional integral in various situa- 
tions, wc have used the properties of linear changes^ of variable of integration 
and properties of functional determinants; they are given in two appendices. 
The transformation S in formula (1.13) is related to the Cartan development 
map; the properties of this map are discussed in the third appendix. 



Linear changes of variable of integration in an infinite-dimensional space 
are sufficiently powerful and varied for the purposes of this paper. In a later 
publication we shall present nonlinear changes, simplifying and generalizing 
earlier works such as [12]. 
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II - A general theorem 



The primary goal of this section is to define functional integrals over L^'^ 
pointed paths (paths with a fixed point) taking their values in a manifold A^, 
by reducing them to functional integrals over paths taking their values in a 
flat space W^. The main definition is given as follows 

(t/T0) (xo) / Vs,Q,z ■ exp (--Qo(^)) <P {xo ■ S(T, z)) . (II.l) 

All the notations are given in paragraph 1. The previous integral, being also 
denoted ^{T,xo), is a solution of the generalized Schrodinger equation 

with initial condition \E'(0, xq) = 4> {xq). 

We give also a general construction of time-ordered products in the form 
of a functional integral generalizing equation (II.l), namely 

([/#(/)) (xo) = / V,,Q^z-e^v(--Qo{z))F{T,z)4>{xo-T.{T,z)). (II.3) 

^ZoJ ^ ^ ^ 

In paragraph 2, we construct the simplest functional integral of type (II.l) and 
prove that it satisfies equation (II. 2). In paragraph 3, we study the general 
case. 

1. The setup, and a summary of notations. 
1.1. A manifold and vector fields. 

- A finite-dimensional manifold N . 

- One-parameter groups acting on A?", denoted cr(a)(r); here a takes the 
values 0, 1, 2, • • • , d and r is a real parameter; the transform of a point x in N 
under (T(Q,)(r) is denoted by x ■ (T(Q,)(r), and 

f^(a) (r) O CT(c,) (s) = CT(c,) (r + s) . (II.4) 

- The generator of {c(a)('')} is the vector field X^q,) in N such that 

d 

— {x ■ (7(«) (r)) = {x ■ (7(«) (r)) (II.5) 
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and in particular 

r=0 

for any point x in N. We do not assume that the vector fields ^(q) commute, 
hence [-^(a); 7^ in general. We often write Y for X(o) emphasizing its 

special role. 

- JCx denotes the Lie derivative w.r.t. the vector field X. 

1.2. Pointed paths on the flat space M.'^. 

- T is a time interval of length T, hence 

T=[ta,h] , T = h-ta. (II.7) 

- to is a chosen element of T; the standard choices are Iq = ta or to — ft- 

- Zq (or Zo,T if we need to specify T) consists of the real vector-valued 
functions z — (^^, • • • , z'^) whose components t 1— > z°'{t) are continuous func- 
tions with square- integrable derivatives i". We assume the normalization 

- Zq is the space dual to Zq. Its elements are interpreted as vector- 
valued distributions z' = {z'l, ■ • • , z'^), each component being the derivative of 
an L^-function. The duality is given by 

{z\z)^ f dtz'^{t)z^{t) (II.8) 

JTJ 

(summation over a and integration over t). 

- s is a parameter equal to 1 or z; its square root ^/s is normalized as 
follows 

f 1 if s = 1 



e'^*/^ if s 



- h = {ha/3) is a constant invertible symmetric real matrix of size d x d. 
We denote by (/i'*^) the inverse matrix and we assume that h is positive 
definite in case s is equal to 1. By a suitable linear change of coordinates, we 
may take h into a diagonal form 

/i = diag(l,---,l,-l,---,-l) (11.10) 

with p elements +1, q elements —1 and p + q — d. Hence p — d, q — in the 
case s = 1. 

- We introduce a quadratic form Qo on Zq as follows 

Qq{z)= I dthapz^{t)zf^{t). (11.11) 

JTJ 
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The corresponding kernel is given by 



DaB\u,r) = -^- — , , ^ ^ , , , (11.12) 

hence the representation 

Qo{z)^ [ du [ drDo,/3{u,r)z"{u)z^{r) . (11.13) 
Jt Jtj 

- On Zq we consider a quadratic form 14^0) with kernel G"^{u,r). As 
above, we have the relations 

c^-vr)=^ ''^f;;. , (11.14) 

2 5^^(w)5^^(r) 

Wo(^')= / / drG'"^(w,r)4(w)4(r). (11.15) 

- We assume that the quadratic forms Qo on Zq and VFq on Zq are inverse 
to each other in the sense of relation (1.6). In terms of kernels, this is expressed 
as follows 

/ dt D^p (si, t) G^^ (t, S2) = SZ S (si - S2) . (11.16) 
Jtj 

Here are explicit formulas for the kernels: 

Dai3{u,r)^ [ dth^f3S'{t-u)6'{t-r)^-h^f3S"{u-r), (11.17) 

that is D : Zq ^ Zq is the differential operator with matrix ^— /ia/S^^ 

Hence, G"'^(w, r) is the corresponding Green's function taking into account 
the boundary condition z'^{to) — 0, that is 

jjOcP (^ _ 7- _ ^g) for u > to,r > to, 

C'f^iu, r)= { /i"^ inf {to -u,to- r) for u<to,r< to, (H.IS) 

otherwise. 



Remark. We shall refrain from integrating by parts in formulas like (H.ll) in 
order not to have to make explicit statements about boundary conditions. 

1.3. Functional integrals on VqW^. 

- The space of paths Zq shall also be denoted by VqW^ to remind us of the 
flat space R*^ where the paths lie and of the fixed point (origin in R^) of the 
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paths z. We shall consider a variety of pairs (Q, W^) consisting of a quadratic 
form Q on a space Z and a quadratic form W on its dual Z' satisfying the 
analogues of relations (11.12), (11.14) and (11.16). 

- For such a pair (Q, VF) we have a translation invariant integrator Vs^qz 
on Z characterized by 

/ Vs,QZ ■ exp (--Q{z) - 2TTi {z', z)) = exp(-7rsW^(z')) (11.19) 
Jzj ^ ^ ' 

for z' in Z'. The normalizations are chosen so that Gaussian integrators and 
Fourier transforms do not include powers of tt depending on the dimension of 
their domain of definition. 

We also write 

Vi^s,Q{z) = V,,QZ ■ exp (-^<5(^)) (11.20) 

and refer to both V^^q and T>liJs,q as "Gaussian integrators". When working 
with the basic pair {Qq, Wq) we simply write Vs and Vujs. In the context of an 
application where s has been chosen once for all, we omit it in the notations. 

1.4. Functional integrals on Vx^N . 

- We fix a point xq in N and consider continuous paths x :T ^ N with 
the fixed point x{to) = xq and square-integrable velocity^. The set of all such 
paths is denoted by VxqN. 

- The time interval T being given, consider an element 2; of Zq. As we 
shall see in Appendix C, the differential equation 

r dx{t) = X^oc){x{t))dz'^{t) + Y{x{t))dt 
\ x{to) = Xo 

admits a unique solution x{-) in VxqN. 

- The previous construction defines a parametrization P of the space 
VxqN of pointed paths on N by the space VqW^ of pointed paths in E*^, 
namely 

P : -p^-^d ^ p^^^ 

by taking z into x. If necessary we shall denote by x{t^ z) the solution of the 
differential equation (11.21) for given z in VqU^. Hence x{t^ z) is a function of 
t and a functional of z. 



^ More precisely, for every smooth function / in C'^{N), we assume that 
the continuous function t f{x{t)) on T is the primitive of a function in 
L2(T). 
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- Assume now that T = [0, T] and to = 0. With the previous definitions, 
define T,{T,z) as the transformation taking a point xq in N into the point 

x{T,z). 

- Take a good^ function (j) on N. Define a functional (p on Zq by 

It is integrable under the integrator Vu>s on Zq. By integrating we get 
/((/), T, xo) = / Vu;s{z)(j){z) 

VsZ ■ exp (^--Qo(-2)) ■ K • S(T, z)) . 

We could replace Qo by another quadratic form on Zq. 

- The functional operator Ut associates to the function (j) on N the 
function xq i— > I{(f), T, xq) on N. 

1.5. Variational techniques. 

A frequently used technique consists in introducing one-parameter vari- 
ations in the space of paths, or in the space of functionals on the space of 
paths. We work in particular with the following variations. 

(i) Fix a time interval T = [t^, ti,] and consider the space VW^ of paths 
x : T ^W^, say of class^ C^. Assuming an action functional S : VM. — M 
(for instance, the time integral of a Lagrangian L{x,x,t)), the critical points 
of S will form a 2(i-dimensional manifold U^^, the so-called space of (classical) 
motions, denoted by Xd- We can parametrize them by a set of parameters 
H = [fi^, ■ ■ ■ , /U^'^); for a given k between 1 and 2d, the derivative dxc\{iJL) /dfJ' 
defines a variation through classical paths, and wc get 2d of them. 

(ii) If to is a given epoch in T, the pointed paths arc defined by the 
condition x(to) = 0. They form the space VqW^. In paragraph III.l we shall 
use a one-parameter family of pointed paths x{X) in (for A running over 
[0,1]). 

(iii) Introducing d independent boundary conditions at ta, and similarly 
at th, we define the subspace Va,b^ of VW^. In paragraph B.2, we consider 
a one-parameter family of paths x{\) in Va,b^i with A in [0, 1], such that 
a;(0) belongs to U^'^ fl Va,b^- The presence and nature of the caustics is 
conveniently analyzed in terms of this intersection. 

(iv) In paragraph III. 2 we consider a one-parameter family of action func- 
tionals S {v) defined on a space VqW^ of pointed paths. 

^ Any function in C°°{N) with compact support will do. 
^ That is, continuous with continuous first-order derivatives. 



(11.22) 
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2. The one-dimensional case. 

We begin by the simple case d = 1. Here X = is a vector field on 
N generating the one-parameter group of transformation (7(r) on N. Hence 
a{r) obeys the differential equation 

d{xo ■ a{r)) = X{xo ■ a{r)) dr (n.24) 

for a fixed xq in A^. We set Y = and consider the time interval T = [0, T] . 
The differential equation (n.21) reads now as 

f dx(t) = X(x(t)) ■ dzit) 

I (n.25) 

x{Q) = Xo . 

Substituting r = z{t) in (n.24) we see that the solution to the previous 
equation is given by x{t) = xq ■ (T{z{t)). Hence the transformation T,{T,z) 
is simply a{z{T)). 

The path space Zq consists of the L^'^ functions z : [0, T] — > R such that 
z{0) = 0. It is endowed with the quadratic form 

Qo{z) = I dtz{tf . (n.26) 
Jo 

The corresponding integrator is, for simplicity, denoted by VgZ, hence our 
basic path integral specializes to 

*(T,a;) = j VsZ- exp ^-^ ^ dtz{t)'^ ^{x ■ (t{z{T))) . (n.27) 

Our goal is to show that we have solved the differential equation 

^^{T,x) = ^C\^{T,x). (n.28) 



Fix a point x in N and define a function of a real variable h{r) = 
(f){x ■ cr{r)). From the group property cr(ri) o a{r2) — cr{ri + r2), we get 

(j){x ■ a{r) ■ a{z{T))) = h{r + z{T)) . (11.29) 

Then define 

H{T, r) / VsZ ■ e-''^°^^^/'h{r + z{T)) . (H.SO) 
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Since the integrand h{r + z{T)) depends on the path z through z{T), a hnear 
change of variable z i-^ z{T) transforms immediately this functional integral 
over Zq into an ordinary integral over M (see formula (A. 38) in Appendix A). 
It is easier to use directly the properties of Fourier transforms of Gaussian 
integrators underpinning (A. 38). Denoting by h the Fourier transform of h, 
we obtain 

h{r + z{T))= [ dph{p)e'''''P^''+<^^K (11.31) 
Therefore, after changing the order of integration, we get 

H(T,r)^ [ dph(p)e^'''P'^ [ V,z ■ e^p(--Qo(z) +2nipz(T)) . (11.32) 

By the definition (11.19) and the duality pz{T) — {pSx, z), we get 

/ VsZ ■ exp(--Qo{z) + 2T:ipz{T)) = exp{-nsWo{p6T)) ■ (11.33) 

By (11.15) and (11.18) we obtain 

Wo{p6t) = G{T, T) = p'T . (11.34) 
Collecting equations (11.32) to (11.34) we conclude 

H(T,r)= f dph{p)e^'''P'-'''P^^ (11.35) 

m 

and by derivation under this integral sign, we conclude 



dH s d^H 



dT 471 



(11.36) 



The vector field X is the generator of the one-parameter group of trans- 
formations cr(r). Hence, for every integer m > 0, we get 

-j h{r + z{T))^{C^cl^){x-a{T + z{T))). (11.37) 

By differentiating under the integral sign in formula (11.30), we deduce 

— j H{T,r)^ j^VsZ-e^p{--^Qo{z))c^<l){x-a{r + z{T))). (11.38) 
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Using definition (11.27), we get therefore 



|:) H{T,r) 



= £^*(T,a;) (11.39) 

r=0 



by applying m times the diff'erential operator Cx acting on functions of x. In 
particular, for m = 0, we get 

HiT,r)l^, = ^iT,x) (11.40) 

and both sides can be derivated with respect to T, giving 



d 

dT ^ ' ' 



= ^^{T,x). (11.41) 



dH s d"^ 

Setting r = in the relation = —-^-^H and using relation (11.39) for 
m = 2, and relation (11.41), we conclude the proof of the differential equation 

A*(T,.) = ^/:^*(r,x). 



Remark. The previous proof can be recast in the following symbolical way 
[5]. We begin with a consequence of equations (11.33) and (11.34), namely 

/ VsZ ■ exp(--Qo{z)) exp{2nipz{T)) = exp(-7rsp^T) , (11.42) 

one of many characterizations of our integrator VsZ. Make the formal substi- 
tution 

and get 

^ VsZ ■ exp(-^Qo(^)) exp(^(T)£x) = e^p(^^Cj,^ . (11.43) 

Apply this operator identity to the function (f){x) and notice that the operator 
exp(r£x) transforms (f){x) into (f){x ■ a{r)). Hence we get 

^ VsZ ■ exp(-^Qo(^)) ■ ct{z{T)) = (^exp ^Cj,^ 0(a;) (11.44) 



that is 



*(T,x) = exp(^^£^ ) 4>{x). (11.45) 
16 



This integrated form is equivalent to the differential equation (11.28) together 
with the initial condition 

*(0,x) (11.46) 



3. The general case. 

3.1. The group property. 

The functional operator Ut defined by formula (11. 1) satisfies the group 
property 

Ut o Ut' = Ut+t' (11.47) 
for T > 0, T' > 0. The proof rests on three facts. 

a) Group property for the point-transformations S(T, z): 

The relevant function space Zo,t consists of the paths z : [0,T] 
with derivative. For z in Zo,t and z' in Zo,t', we define a new path z x z' 
by the rule 

f z(t) for < t < T 

(z X z'm = 1 (11.48) 
(z{T) + z'{t-T) for T<t<T + T'. 

It is obvious that 2; x 2;' is an element of Zo,t+T'- Furthermore, the map 
(2;, 2;') I— > 2; X z' is an isomorphism of Zo,t x Zo,t' onto Zq, t+t' and, by the 
uniqueness of the solution of the differential equation (11.21), we obtain 

xq ■ E{T + T',zx z') = xo ■ S(T, z) ■ S(T', z') (11.49) 

(see also formula (C.23) in Appendix C). 

b) Quadratic forms: 

We denote by Qo,T the basic quadratic form on Zo,t given by equation 
(11.11), namely: 

Qo,t{z)= I dtha,i3z'^{t)z^{t). (11.50) 

Using similar definitions for (5o,t' and (5o,t+t', one obtains immediately 

Q0,T+T'{Z X Z') = Qq,t{z) + Qo,T'{z') (11.51) 

and by exponentiating 
exp(^-^(5o,T+T'(^ X z')^ = exp(^-^(5o,r(^)) exp(^-^(5o,T'(^')) ■ (11.52) 

We can identify Zo,t+t' to Zo,t x Zo,t'- This identification entails an 
identification of the dual space Zq t+t' product space Zq ^ x Zq ^, . 
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We use the notations ( for elements of Zo,t, C' in Zo,t' and denote hy ( x (' 
the corresponding element in Zq j'^rp, . Using simple algebra, one derives the 
identity 

Wo,T+T'{C X C') = Wo,t{0 + Wo,T'{C) (n.53) 

from (11.51). Here Wo,t denotes the quadratic form on Zq ^ inverse to Qo,t, 
etc. 

c) Integrators: 

The following form of Fubini's theorem holds 

f VsU ■ F{u) = /" VsZ I Vsz' ■ F{z X z') . (11.54) 

According to the general method explained in Appendix A, it is enough to 
check this formula for a function of the form 

F{z X z') = exp(-27ri((C, z) + (C', z'))) . (11.55) 

But then our contention follows from the characterization (11.19) of the inte- 
grator, and the relations (11.51) and (11.53). 

We combine this formula with equation (11.52) and obtain another form 
of Fubini's theorem 

f Vus{u)G{u)= I Vuj,{z) ! Vujs{z')G{zx z'). (11.56) 

We conclude the proof of equation (11.47). Indeed 
{Ut{Ut'<P)) {xo) = [ Vus{z) {Ut'<P) {xo ■ S(T, z)) 

= 1 Vuj,{z) I Vu;s{z')(j){xo-^{T,z)-J:{T\z')) 
= / VuJsiz) / Vuj,{z')(j){xo-J:{T + T',zx z')) 
= / Vco,iu)(j){xo-^{T + T',u)) 

= {UT+T'(j)) {xq) . 
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3.2. The differential equation. 

^From the group property (11.47) it follows that we need to establish the 
partial differential equation (II. 2) at the time T = 0, and that the general 
case will follow. That is, we want to prove 

(C/t0) (xo) = 0(a;o)+T(^/i"^£x(„)/:X(,)</'(a;o) + >Cr0(a;o)) +o{T) . (11.57) 
This relation implies also the initial condition 

hm{UT4>){xo) = <t>{xo) . (11.58) 

For the proof, we rely on the scaling properties of paths as described in 
paragraph A. 3. 6. For z in Zo,i we define the scaled path zt in Zo,t by 

ZT{t) = T^'^z{t/T) . (11.59) 

By scaling the differential equation (11.21), we get 

dx{t/T) = T^f^X^^){x{t/T))dz^{t) + TY{x{t/T))d{t/T) . (11.60) 

Hence the transformation S(T, 2;^) in N defined using the vector fields X(q,) 
and Y is the same as the transformation S(l, 2) using the vector fields T^/'^X(^a) 
and TY . Moreover we can transfer the path integration from the variable 
space Zo,T to the fixed space Zo,i, that is 

{Ut<I>){xq) = j Vu;s{z)(f>{xo ■ ^{T,zt)) . (11.61) 

We use then a limited expansion of (f) around ^(xo), namely 

4>{xo ■ E(T, zt)) = 4>{xo) + T^/^Cx,^^cl){xo)z"{l) 

+ |£x(„)>Cx(,)0(xo)^"(l)^^(l) (11.62) 
+ T£y0(a;o) + o(T3/2) . 

Recall the integration formulas 

Vu;s{z)z'^{l) = (11.63) 



0,1 



and 

Vu;s{z)z''{l)z'^{l) = s/t"^/27r. (11.64) 



^0,1 
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Collecting equations (11.61) to (11.64), we conclude the proof of (11.57). 

The previous calculation can be extended to give the complete Taylor 
expansion of (?7t0) (xq) around T = 0. Since the vector fields Xf^^) ^nd Y 
do not commute, we have to use a time-ordered exponential to express the 
solution of the differential equation (11.21). 

3.3. About the construction of S(T, 2;). 

For simplicity, assume Y = and hc^p = Saj3- Consider again a path 
z : [0,T] of class L^'i with ^(0) = 0. For each t in [0,T], denote by 

z^{t), • • • , z'^{t) the components of the vector z{t) and define the transforma- 
tion s{z{t)) on N by 

xo ■ s{z{t)) = xo ■ a^i){z\t)) • . . . • a^d){z'^{t)) . (11.65) 

Here {<7(a)(^)} denotes the one-parameter group of transformations in N with 
generator X(^^y In general, the vector fields -'^(a) do not commute, hence the 
transformations s{z{t)) do not form a group. In the general case, we can use 
a multistep method to solve the differential equation (11.21), hence 

xo-E(T,z) = lim xo-s{z(T/n))-s{z{2T/n)-z{T/n))-. . .•s(^(r)-^(r-T/n)) . 

n=oo 

(11.66) 

Putting this into the integral (11. 1), we obtain after some calculations the 
following variant of the Lie, Trotter, Kato, Nelson formula 

^^"^ =iL^ ■ ■ ■ 4%) ^ <!> ■ (n-67) 

Here C/^"^ for a in {1, ■ ■ ■ , d} corresponds to the one-dimensional case studied 
in paragraph II. 2, hence 

C'^"'=e-p(i£i<.,)- (11.68) 

These relations are in agreement with 

4. Some generalizations. 
4.1. Including a potential. 

Consider the Schrodinger equation with potential 

^ = 7z/^"'^^x,.,/:x,,,* + F*. (11.70) 
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A path integral solution is obtained as follows 



^{T,xo) = j VsZ-e^^{^'^Qo{z) + J dtV{xo ■ ^{t, z))^ 4>{xo ■ ^{T, z)) . 

(11.71) 

The proof can be obtained by a slight generalization of the arguments pre- 
sented in paragraph II. 3. We can also use the following trick: we add one 
variable ©, considering the manifold x M. The vector fields A^^,) in N give 
vector fields, also denoted by A^, in A/" x M, which have a zero component 
on the factor R. Moreover Y = V{x)d/dQ. The function *(T,a;) satisfies 
the equation (11.70) if, and only if, the function \I/(T, x) exp O satisfies the 
equation (II. 2) on A" x M. The differential system (11.21) is now written as 



dx — X(^f^-^{x)dz'^ 
d@ = V{x)dt . 

Hence the transformation S(T, z) takes (a;o,0o) into 



(11.72) 



a;o-S(T,^), Go + / dt V{xq -^.{1, z))] , 



T 



hence 



^{{xo,Q) ■T.{T,z)) = (t){xo-T.{T,z))e^v{^j^dtV{xQ-T.{t,z))j . (11.73) 

Equation (11.71) follows easily from these remarks. Notice that the exponent 
in this formula does not have a simple dynamical interpretation in general - 
but see our illustrations in section IV. 

4.2. Time-ordered product. 

Suppose that F is a suitable functional on the space Zo,t- We generalize 
equation (II. 1) as follows 

{U^cf) {xo) = I Vusiz) F{z) (j){xQ ■ E(T, z)) . (11.74) 

By imitating the proof of (11.47), we get 

U^^T, = U^U^', (11.75) 

for functionals F on Zo,t and F' on Zq^t', where the functional G on Zo,t-i-t' 
is defined by 

G{z X z') = F{z)F'{z') . (11.76) 
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Ill - Semiclassical Expansions 



We shall compute the semiclassical approximation of the expression in the 
basic equation (11. 1). The scalar potential V is included via coupled equations 
rather than via an additional variable, so that the equations will take readily 
their familiar forms; i.e. we compute the wave function: 



For ready use of our results in quantum physics, we compute Xh)- There- 
fore in the general set up, we set to — ti,, xq — x^, hence Z is the space Z{, 
of paths vanishing at t^ and x{ti„ z) = Xb- In paragraph 3, we give the corre- 
sponding results for '^{tb,Xa)- 

We choose an initial wave function given by 



where T is a smooth function on N with compact support, and Sq is an 
arbitrary, but reasonable function on N. 

The initial wave function given by (III. 2) generalizes plane waves on 
W^, but is, obviously, not a momentum cigcnstate. We can nevertheless 
call the semiclassical expansion of \l/(t5,X5) with this initial wave function 
a momentum-to-position transition amplitude" for three reasons: 

(i) In the limit h = 0, assuming s = i, the current density corresponding 
to the initial wave function 4> is 



Consequently, \I/(t&, Xb) is the amplitude corresponding to the transition from 
momentum VSo{x) to position Xb- 

(ii) We shall expand (IH-l) around a "classical path" z^i characterized 
by initial momentum and final position. 




(III.l) 




(111.2) 
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(iii) The leading terms (henceforth labeled WKB) of semiclassical ap- 
proximations combine as if they were transitions between eigenstates, e.g. 

(position I position)^j^g ~ J momentum (position | momentum^^^^^ 

(momentum | position)^j^g 

[see reference 15, §7 in the Appendix]. 

Another useful initial wave function is 

0(x) = 4„(^). (111.4) 

The corresponding expression '0(tb, xi,) gives a position-to-position transition 
amplitude. In general, this case is more complicated than the previous one, 
because now the initial wave function restricts the domain of integration Z. 
Moreover the paths z E Z such that x{ta, z) — Xa do not usually have a 
common origin. Even at the very best, when N = W^, and z (ta) ^ Xa — Xb, 
one needs to be careful because the domain of integration, say Za,b C Z, is not 
a vector space but an affine space. Therefore, we shall compute position-to- 
position transitions in section IV where we specialize the basic equation (III.l). 
We shall treat in detail the case N = and give the necessary indications 
and references when N = 0(M), a frame bundle over a Riemannian space 

Semiclassical expansions are best analyzed in the broader context of 
spaces of paths with no requirement on their boundary values. For instance 
let VM.'^ be a space of paths z with no requirement on z (ta), nor z (t^), and let 
Va,b^'^ be the subspace of VM!^ with d requirements at ta and d requirements 
at tb- Let {S{h',z)}j^ be a one-parameter family of actions, and let U'^'^{v) 
be the 2(i-dimensional space of motions made of the critical points of S{v, z). 
A classical paths Zc\ G Va,b^ for the action S{y^ z) is at the intersection of 
Va,b^ with U'^'^^f). If the intersection is transversal, no Jacobi field of Zc\ is 
in the tangent space to Va,b'^- In this paper we assume that such is the case. 
Otherwise there are caustics and we refer the reader to the literature [e.g. 2] . 

1. General strategy. 

We introduce a Lagrangian 

L(t, z) = ^haf3Z^it)zl^{t) - Vixit, z)) (III.5) 

where as usual x{t,z) is the solution of the differential equation (11.21) with 
boundary condition x{tb, z) — Xb- It is a function of z{t) and a functional of 
z. From this Lagrangian we deduce the action functional 

S{z)^ I dtL{t,z)-\-So{x{ta,z)) (III.6) 
Jtj 
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on the space = Vo^'^ of paths z obeying the boundary condition z{tb) = 0. 
Let Zci in be a critical point of the action functional S. 

As it is customary in the calculus of variations, we take a one-parameter 
variation 

^(A) = ^ci + AC (III.7) 

with C in Zb and the equation 

= (IIL8) 



^5(.(A)) 



A=0 



has to be satisfied for all (. That yields, after integrating by parts, a functional 
differential equation for Zc\{ta)- 

Under the affine change of variable from z to ( given by (III. 7), we obtain 

Qoiz) = Qoiza + AC) = Qo(^ci) + 2AQo(^ci, C) + >^^QoiC) (ni.9) 

and 

Vs,Q,z = P,,A^Q„C- (ni.io) 

The expansion of x{-, z{X)) around x{-, Zc{) reads 

P (zci + AC) = P (^ci) + AP' (zci) • C + ^A^P" {z,i) ■ CC + 0(A3). (III.ll) 

Here P' (zd) and P" (zd) are the first and the second derivative mappings of 
P at Zc\, where P : Z^ — > Vx^^ takes z into a;(-, z). We abbreviate x{t^ Zd) to 
Xc\{t)- They are of the form 



Jt OZ'^-y{S) 



= - dsk^f3{t,s)C^{s) =: C^(t) 



" ds [ ' duk'^i3^{t,s,u)C'^{s)C{u). (III.13) 
t Jt 



We could obtain the short time propagator by expanding the quantity 
(f) (xb • z)) in equation (III.l) around cf) (xj,) with x^ e N. Here, we shall 
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expand (xb ■ z)) around (f) {xb ■ Zd)) with Zc\ G Zb. When we choose 
for Zci a critical point of the action, we obtain the WKB approximation. 
These are two different expansions, and in general, the short time propagator 
is different from the WKB approximation [14] . But in some simple cases, they 
are equal [1]. 

Remark. A simple and powerful argument of Stephen A. Fulling clarifies 
this issue: The Schrodinger operator exp{—itH/h) can be written in terms 
of dimensionless terms ^ = -^AA + BV, where A = ^ and B = ^ with 
A a coupling constant. There are 4 different possible expansions of physical 
interest: 

- expansion in B, equivalently expansion in A; 

- expansion in A, equivalently expansion in m~^; 

- expansion in A - B, equivalently expansion in t; 

- expansion in A/B, equivalently expansion in h. 

2. Momentum-to-position transitions. 

In this paragraph we take s = i and write I^s,go simply as ^^Qq. We 
expand the integrand in (III.l) by taking z in the form Zd + C where Zd is a 
critical point of the action functional S, with boundary condition Zd (tb) = 0. 
The initial wave function is given by (III. 2). 

The terms independent of ( combine to make the action function 

S{tb,Xb) = ^Qo{zd)- [ ' dtV{xdit))+Soixdita)) (III.14) 
where Xdit) is by definition x (t, Zd-)- The terms linear in ( 

dt (haf3Z^,{t)z^,{t) - VaV {xd{t)) {P' (zd) ■ (t)) 



L 



(III.15) 



+ Vo.So{Xd{ta)){P' {Zd)<T{ta) 



vanish since Zd is a critical point of S such that Zd{tb) = 0. The quadratic 
terms can be written in the form {Qo + Q) {(,)■ Thanks to the equations 
(A. 127) and (A. 132), the corresponding integration is: 

^PQ„C-exp(-^(Qo + Q)(C)) =Det(Qo/(Qo + Q))'/' • (III.16) 

Thus the dominating term of the semiclassical expansion of '^{h, Xb) is: 

*WKB (tb, Xb) = exp (^-^Sitb, Xb)^ ■ Det (Qo/(Qo + Q)f''^ ■ T {xd{ta)) ■ 

(III.17) 
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We proceed to calculate Det{Qo/Qv'), with Qi, = Qo + ^Q, in the case 
V = 1. Physical arguments as well as previous calculations (see [6] and the 
appendix of [15]) suggest that Det((5o/Qi) is equal to detK^^j where 
the Jacobi matrix K'^ p {ti,, ta) is defined by: 

K%{t,^t.)^^^^ (111.18) 

Therefore we shall construct a one-parameter family of actions S{v, z) and a 
one-parameter family of Jacobian matrices: 

K",(.;4,t«):=|fj^. (III.19) 

[v; ta) 



We shall show that: 

det (^K{iy; ta) ■ K{0; t^, ta)'^) =: c{iy) (III.20) 

satisfies the same boundary condition as Det {Qq/Qv)i namely 

c(0) = Det {Qo/Qo) = 1, (111.21) 
and the same differential equation as Det [Qo/Qv)-, namely, 
d 



, -in 
dv 



Det (Q./Qo) = (^Q;1 Aq,^ = Ty{Q-^Q) (III.22) 



(See equation (A. 124) in Appendix A). 

2.1. Differential equation satisfied by c{i'). 

We choose the family S{i', z) such that its second variation evaluated at 
Zci be 

s" {u, z,0 ■ cc = Qo(C) + i^QiO QAO ■ (ni.23) 

It follows from its definition (III. 18) that the /^-column Kfj [v^t^ta) of the 
Jacobi matrix K{v^ t, ta) is the following Jacobi field along the classical path 
Zciiv-i t) of the system governed by the action S{v^ z) 

Q,K^^^{u-t,ta)^0 (III.24) 

K^^f3^{iy;ta,ta) = S^. (III.25) 
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Relation (III. 24) encodes a functional differential equation for K and an im- 
plicit equation for dK{y\ t, to) /dt aX t — ta- It also gives, after differentiating 
w.r.t. v 

Q^-^K{u) = -QK{u) (111.26) 
which can be solved with the retarded Green's function of Qi,, namely: 

^K{iy; t, ta) = -J^ GT{t, s)Q{s)K{v- s, ta) ds . (III.27) 

The retarded Green's function can be expressed in terms of Jacobi fields, 
namely: 

Gl^\t, s) = 9{t - s)J{v; t, s) (III.28) 

where 6 is the step function, equal to 1 for t > s and for t < s; more- 
over J^fjtjs) is a Jacobi matrix, i.e. each column is a Jacobi field. The 
a-component of the /3-Jacobi field J^^^ {v; t, ta) is 

J-^'H-;t,ta) = (III.29) 

where Pci{i^,t) := SS{i', z)/Szci{i',t), here Pci,a{i^,t) = hai3Z^^{i>,t). 
^From (III. 20), we get the equation for c{v): 



d 

dv 



Jin c{v) = tr (^N{i^; ta, tb) -^JK{i^; tb, ta)^ (III.30) 



where the matrix N is the inverse of K 

N{iy;ta,tb)K{iy;tb,ta) = l. (III.31) 

Using (III. 27) and (III. 28), we obtain the sought-for differential equation: 
-^jenc{u) = -ti(^J dsJ{u;tb,s)Q{s)K{u;s,ta)N{u-ta,tb)^. (III.32) 

Within focal distance of Zc\ {i^',ta), the d Jacobi fields K(^)(i/) and the d 
Jacobi fields j'^^\i>) are linearly independent, and c{v) is defined. It gives the 
rate at which the fiow of classical paths {zdiy)} diverges or converges. 
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2.2. Comparison between the two differential equations (III. 32) and (III. 22). 

We shall express in terms of Jacobi fields, as this is always possible 
within focal distance of Zc\ (to)- The unique inverse G^, of Qi, satisfies the 
equation 

QuG^ = t (III.33) 

with 

Gy{tb,s) = Q (III.34) 

by virtue of (A. 40) together with the specialization used in (A. 62) when ap- 
plied to the space Z;, of paths vanishing at ti,. Therefore (see formula (B.26) 
in Appendix B) 

G,{t,s) = d{s-t)K{v-t,ta)N{v-ta,tl)J{v-tk,s) 

^ ^ (III. 35) 

- e{t - s)J{u; t, tb) N{u; t^, t^) K{v- ta, s) 

where K is the transpose of K, and N is K's inverse. 
Substituting Q'^ = G^ in (III.22) gives: 



^ J£n Det(Q^/Qo) =tr( j dt K{iy;t,ta) N{iy;ta,tb) J{iy;tb,t) Q{t) 



du 



(III.36) 



Comparison with (III. 32) shows that 

£n Det {Qu/Qo) = -£n c{u) (III.37) 



i.e.io 



Det (Qo/Q^) = det (^K(u; tb, ta) ■ K(0; tb, ta) ') • (111-38) 
Here K{0; t, ta) satisfies 

QoK{0;t,ta) = (III.39) 
K{0;ta,ta)^l (III.40) 

where 

Qo{z) = -j^dt haf3Z^{t)zf^{t) - haf3Z^{ta) {ta) • 



The credit for the technique used in deriving (III. 38) is due to B. Nelson 
and B. Sheeks [6]. It is made simpler and more general here by not integrating 
Qi, by parts. The first calculation giving the ratio of functional determinants 
in terms of finite determinants can be found in [16, 5]. 
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Hence we have 



^K{0;t,ta) = , j^K{0;t,ta) 



= (III.41) 



therefore 

K"^(0;t,t„)=5«^. (III.42) 
Inserting Det((5o/Qi/) = detK'*^ (tbita) into (III. 17) we obtain: 

det ^^7^ ) ■exp(^S{h,x,)]-T{x,,{ta)) • (III.43) 

If Zc\ (tb) is conjugate to Zd (ta), in the sense of caustic theory, one needs [2] 
to include terms in of order higher than 2 in the calculation of 

2.3. End of calculation. 

Equation (III. 43) is not the end of the calculation; ^'wkb rnust be ex- 
pressed in terms of Xci : T — > N, where 

Xci{t) = x{t, Zci) = • S (t, Zci) (III.44) 

but not in terms of : T — > M*^. Two techniques present themselves: we know 
the evolution equation satisfied by the wave function \1/ on N, therefore we 
can construct an action on and find its critical points. But we do not need 
to find Xc\ corresponding to Zc\, we only need the expression corresponding to 
the prefactor of (III. 43). This can be obtained by a simpler technique. 

Recall the parametrization P : VqW^ T^x^N given by P{z) — x. Corre- 
sponding to the action functional S on Vq^^ we get a functional S on Vx^N 
such that S = S o P, i.e. S{x) = S{z). Making the substitution z = + C 
and expanding up to second order terms in we get 

S{x,i) = S{z,i) (III.45) 
S'{x,,) ■ = S'{z,,) (III.46) 

s"{xd)-p^P^ = S"izd). (III.47) 

In the beginning of this paragraph we calculated the expansion of <S'(2;ci -|- C) 
up to second order and obtained the expression 

s"{z,,) ■ cc = Qo(0 + QiO = QiiO (111.48) 
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for the second order variation. We also know that the infinite-dimensional 
determinant Det{Qo/Qi) is equal to the determinant of the Jacobi matrix 

dzci ih) /dz^x {to). 

Similarly we write 

s"{x,{) . ^C^C = Qo(C) + QiO (ni.49) 

and identify Qq as follows. Let 

So{z) := \Qo{z) (III.50) 

and 

So{x{z)) -.^ So{z) (III.51) 

then 

go(C):=:5;'(xe,)-^C^C (111.52) 

det ll^l^ = Det(Qo/(Qo + 0))^ Det(Qo/(Qo + Q)) ■ (111.53) 
Therefore 

det 1^ = det 1^ . (III.54) 

dXc\ [ta) OZcl {ta) 

In conclusion, we obtain the sought-for semiclassical expansion 

*WKB(t6,:Cb) =det(^|^^^^ ^ ■ exp(^^S{tb,Xb)yr{xci{ta)). (III.55) 

According to equation (III. 46), the path Xc\ in A'' is a critical point for the 
action functional S : Vx^N R. By construction, we have 

S\x{-,z)) = S'{z). (III.56) 

Using equation (III. 5) and (III. 6) this can be made more explicit as follows 

S{x)^ \ I dthaf3Z''{t)z^{t)- [ V{x{t))dt + So{x{ta)) . (III.57) 

Jt Jt 

In our general setup , the first integral remains somewhat implicit, but can 
be used for practical calculations in the applications given in section IV. 

The prefactor in (III. 55) also gives the volume expansion or contraction of 
a congruence of classical paths originating in the neighborhood of So{xc\{ta)) 
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with momentum ViSo(m), m E N. This determinant depends both on the 
choice of the initial wave function and the dynamics of the system. 



then 



Set doja the volume element on N at Xci{ta) and dw^ = det ( ] duja', 

lim / \-^{tb,Xb)\^ dujb ^ / \(l){Xa)\'^ dwa 

where Ct belongs to the group of transformations generated by the classical 
flow (see details in [5 p. 299]). 

3. Diffusion problems. 

In a diffusion problem, one is interested in computing the functional 
integral 



^{th,Xa)= J r>s,Qo^exp(^-^Qo(^)) 



1 N (III-58) 



exp( — / dtV{Xa-^{t,z))\ -(PiXa-^itb^z)) 

where Z denotes now the space of paths z vanishing aX t — ta- If we choose (f) 
to be of the form (III. 2), namely 

{xa ■ S(t6, z)) = {x{h, z)) = exp(^-^So {x{h, z))^ T {x{h, z)) (III.59) 

then computing (III. 58), with s = i, can be said to be computing the prob- 
ability amplitude of a transition from a position Xa to a momentum of the 
form 

p (xb) = V5o (xb) (III.60) 
for the end point xj, = x{tb, z) of the path x. 

With (j) given by (III. 59), one obviously expects the WKB approximation 
of (III.58) to be 



5<i (ta) 



1/2 



^'WKB {h, Xa) = det I ^ g ; " exp -S{tb, Xa) ) • T(Xcl(4)) . (III.61) 

Nevertheless, it is gratifying to derive (III. 61) by the method followed in 
paragraph 2. The only necessary changes are described as follows: 

c{u) = det(^K{u;ta,tb) ■ K{0;ta,tby^^ (111.20^^^^) 
Gf^it, s) = e{t - s)J{v- s, t) = -e{s - t)J{iy; t, s) (111.28*'^') 

G^{t, S) = e{s - t)J{v; t, ta) N{V- ta, tb) K{l^; tb, S) g^bis^ 

-e{t-s)K{i^;t,tb)N{u;tb,ta)J{u;ta,s) 
(see formula (B.27) in Appendix B). 
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IV - Illustrations 



The specializations presented in this section illustrate and develop the 
general formulas derived in sections II and III. They differ by the choice of 
the manifold N and of the initial wave function 0. We treat the case of scalar 
wave functions, but the case of a tensor field is easily accommodated by using 
the Lie derivative of tensor fields. 

We are basically interested in the Schrodinger equation. So the reader 
should substitute z to s in the formulas. Moreover, we want to evaluate the 
value of the wave function ^' at the final time ti,, and final position xi,. This 
dictates the choice of the space of paths, determined by z{ti,) = 0. 

1. Point-to-point transitions in a fiat space. 

The basic manifold is the euclidean space of dimension (i, in carte- 
sian coordinates (see paragraph IV.2 for other coordinates). For a vector x, 

with coordinates x^, • • • , x'^ , we set its length to be |a;| = ( X^Q,=i(a;")^ ) as 

usual. We consider a particle of mass m moving in the field of a potential V . 
The Lagrangian and action are given as usual by 

L{x,x) = '^\xf -V{x\ (IV.l) 
5(a;(-)) = y dtL{x{t),x{t)). (IV.2) 

In classical mechanics, we solve the equations of motion with suitable 
boundary conditions: 

mxci = -VV{xci) (IV.3) 

Xcl{ta) = Xa , Xcl{tb) = Xb , (IV.4) 

where Xa,Xb are points in W^. In quantum mechanics, we want to solve the 
Schrodinger equation 

ih— = A* + (IV.5) 

dt 2m ^ ^ 

^{ta,x) = (l>{x). (IV.6) 



with initial condition 



The point-to-point transition amplitudes are given by 

{tb,Xb I ta,Xa) = ^(tb.Xb) , (IV.7) 
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where (f){x) is a delta function S{x — Xa)- 

The original claim of Feynman was that we can solve the Schrodinger 
equation by the following path integral 

^ih, Xb) = [ Vx- e^^(-)/'^ • 4>ixita)) , (IV.8) 

JVb 

where Vb is the space of all paths a; : T — > R*^ with endpoint at Xb, namely 
x{tb) = Xb- The questionable part was the rigorous definition of the integrator 
Vx. 

To fit within our general framework, we consider translations acting on 
M.^, namely 



X 



<^{a){r) = {x^,---,x'^-\x'^ + Xr,x'^+^,---,x'^) . (IV.9) 



The corresponding Lie derivative is given by 

^^(.)/ = ^^ (IV-10) 

and the parameter A is chosen equal to (h/m)^. The general differential 
equation dx = X^^^{x) ■ dz*^ reduces to dx*^ = \dz°^ (for a in 
Hence the solution 

x{t, z) = Xb + Xz{t) (IV.ll) 

describes the parametrization of the space Vb of paths x with x{tb) = Xb by 
the space Z5 of paths z with z{tb) = 0. Accordingly, we obtain for the action 

-\s{x{;z)) = -- f dt\z{t)f + \ I dtV{x{t,z)). (IV.12) 
sn s 7t Jt 

Substituting z for the integration variable x in equation (IV.8), we obtain the 
path integral (remember that s = z in quantum mechanics) 

^itb,Xb) = Vz-exp(^-^S{x{;z))y4>{x{ta,z)). (IV.13) 

^From our general results in section II, the function is a solution of the 
differential equation 

^='yc% * + (IV.14) 

a 
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that is 

dt 2m ^ ^ 

For s = i, this is the sought-for Schrodinger equation. The integrator T>z in 
the space is invariant under translations and is normahzed by^^ 

J Vz-exp(^~ J dt\z{t)f^ =1. (IV.16) 



We derive now the semidassical approximation^^ . We use the initial wave 
function (f){x) — S {x — Xa) and reparametrize the paths around the classical 
path Xci, that is 

x{t,C) = xa{t) + Xat), (IV.17) 

with ( in Zft. Since the integrator in is invariant under translations, and 
Xciita) = we transform equation (IV. 13) into 

^itb,Xb) = j^T^C- exp(^-^5(a;ei + AC)) KK(ta)) ■ (IV.18) 

In expanding S{xci + XQ in powers of A, there is no term linear in A, since 
Xci is a critical point of the action functional S and because C{ta) = Citb) = 0. 
Hence we obtain 

-^S{xa + AC) = -^-S(xci) - - [ dt |C(t)|2 (IV.19) 
sn sn s Jt 

+ — / dtV^Vf3V{xa{t))C{t)C^{t) +0{h'/^) . 

The action S{xci) corresponding to the classical path Xc\ with endpoints 
Xcl{t a) Xq\ {tb) = Xb is nothing else than the classical action function 

S(tb, Xfj, ta, Xa)- 



The scaling factor A = (/i/m)^/^ has been chosen in such a way that no 
physical constant enters in this normalization. With the notations of para- 
graph A. 3. 6, we have the dimensional equations 

[^«] = rV2 , [t] = T. 

Since A = (/i/m)^/^, the semidassical expansion will proceed according 
to the powers of h^^^. 
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1/2 

Omitting terms of order % ' , wc obtain the WKB approximation 
^WKB{tbi xi,) to Using formulas (IV.18) and (IV.19), we derive 

'^WKB{tb,Xb) ^ ex.p(^-^S{tb,Xb;ta,Xa)^ I, (IV.20) 
with the integral 

1 = J^'^C- exp(-^gi(C)) 5{Xata)) . (IV.21) 
Besides the quadratic form 

QoiO^J^dt (IV.22) 
corresponding to the free particle, we need the quadratic form 

Qv{C) = -- I dtKp{t)C{t)C''{t) (IV.23) 

with haf3{t) = 'Va^ f3V{xci{t)). In (IV.21), we use the quadratic form Qi = 
Qo + Qv- 

The easiest method to calculate the functional integral / consists of the 
following steps. 

a) Changing the integrator, according to formulas (A. 126), (A. 132) and 
(A. 134) in Appendix A, we obtain / = /1/2 where 

h = \Bet{Qo/Qi)\'^' ,-ind(QO ^ (IV.24) 
l2= f P.,Q,C-exp(--gi(C))5(AC(ta)) • (IV.25) 

Here Ind((5i) is the number of negative directions for the quadratic form Qi. 

b) Restricting the domain of integration: to treat the S factor in I2, we 
use the linear change of variables ( ^-^ C(^a) from to W^. By a method 
similar to the one used in paragraph A. 3. 8, we obtain 

h = (detsM/l(5^J)-'/'A-^ (IV.26) 

where Wi is the quadratic form on inverse to Qi. Using the Green's 
function G given by equation (B.26), we evaluate the matrix Wi{5t^) as follows 

WiidtJ - G{ta, to) = iV(ta, h) J{h, to) . (IV.27) 
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c) Reducing a functional determinant to a finite determinant: using the 
same strategy as in paragraph III. 2 we obtain 



Det(Qo/Qi) = det (k{1- t^, to) ' K{Q- 1^, . (IV.28) 

Here K{i';t,ta) is the Jacobi field defined by 

Q,K{u;t,ta) = , K{u;ta,ta) = l (IV.29) 

for u equal to or 1. 

Finally, collecting the previous equations (IV. 20) to (IV.29) and using 
equations (B.ll), (B.12), (B.7), (B.8), (B.19) and (III.42), we obtain the 
WKB approximation to the point-to-point transition amplitudeJ (in the case 
s = i): 

{tb,xb I ta,xa)^j,^ = ch-'"^ \deid''S/dxZ ■ dxl^'\'^"' . (IV.30) 

Here S = S{ti,,Xb;ta,Xa) is the classical action function and the phase factor 
c is e'^*'^^"'^-'/^ where p (q) is the number of positive (negative) eigenvalues of 

the van Vleck-Morette matrix (d^S/dx^ ■ dx'^j 

l</3<d 

For the higher-order terms in the semiclassical approximation, whether 
Xci is or is not a degenerate critical point of the action, we refer the reader to 
the literature [5] or the references at the end of Appendix B. The reader can 
easily transfer these old results using the simpler and more general formalism 
of this paper. 



2. Polar coordinates 
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Our basic integral is formulated in terms of a transformation E(t, z) on 
a manifold N, in a form valid for an arbitrary system of coordinates. Before 
considering the case of a general Riemannian manifold in paragraph IV. 3, we 
consider polar coordinates in a plane. The case of cylindrical coordinates in 
is very similar. 

"Path integrals in polar coordinates", the 1964 paper by S.F. Edwards 
and Y.V. Gulyaev [17] has been, and still is, at the origin of many investi- 
gations. To the best of our knowledge, all the papers on the subject deal 
with the discretized version of the path integral, and propose various path 
integral prescriptions when the discretized values of the paths are expressed 
in coordinates other than cartesian. 



Contributed by John La Chapelle. 
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The basic manifold is N = M^\{0} with coordinates and x^. Wc also 
consider the manifold N =]0, +oo[xM with coordinates r,9 (sole restriction 
r > 0) and the covering map II : N —>■ N taking (r, 9) into (^x^, x^) with 

a;^ = rcos^ , a;^ = rsin^. (IV.31) 

Two points of map onto the same point of N if, and only if, their 6- 
coordinates differ by an integral multiple of 27r. 

Let A = (0%)'^ (8%)^ Laplacian in cartesian coordinates. We 

know how to solve the Schrodinger equation 

s 

V = i^A* (IV.32) 
dt A-K ^ ^ 

by means of the path integral 

^{h,Xb)= [ P^-e-'^^°("^/Xa;(ta,^)) (IV.33) 

JZb 

provided that x{t, z) is a solution of the differential system 

dx^ = dz^ , dx^ = dz^ (IV.34) 

with boundary condition x{tb, z) = xi,. To solve the Schrodinger equation 
in polar coordinates, we need only to transform the system (IV.34) in polar 
coordinates with the help of the transformation equations (IV.31). Hence 

dr = cos 6* • dz^ + sin 6* • dz'^ = X^^^dz^ + X^^^dz'^ 

■ ft ft (I^-^^) 

de=-^.dz^ + ^-dz^ = X^^dz^ + Xl^^dz^ . 

The vector fields and X^2) can be read off from the above equations: 

a sin ^ a m7Qft^ 

r d cosO d 

Hence, if (r(t, z), 6{t, z)) is the solution of the differential system (IV.35) such 
that r(t5, z) = Tb, 9{ti„ z) = $1, the path integral 

n, Ob) / Vz- e--Q°(^)/Xr(t«, z),e{ta, z)) (IV.38) 
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solves the Schrodinger equation (IV. 32). As expected, the operator 

is the Laplacian on ]R^\{0} in polar coordinates. 

We want to evaluate the point-to-point transition amplitudes in polar 
coordinates, denoted by (t6,r6,^b | ta,ra,Oa)- To obtain them, it suffices to 
put 0(r, 9) = S{r-ra)S{9-9a) in equation (IV.38). Solving the system (IV.35) 
is easy by reverting to cartesian coordinates, hence 

r(t, z) cos 9(t, z) ■= Th cos 9h + (t) 

(IV.39) 

r{t, z) sin 6'(t, z) — ri, sin 6*^ + z^ {t) . 

For given x^, x'^, the equations (IV. 31) in r,9 have infinitely many solutions, 
and we derive easily 

S{r cos 9 — Ta cos 9a) 5(r sin 9 — Va sin 9a) = — ^ S{r — r^) S{9 — 9a — 2nn) . 

^« neZ 

(IV.40) 

Substituting r{t,z) to r and 9{t,z) to ^ and taking into account equations 
(IV.39), we obtain the path integral representation 

— ^ {tb, n, 9b I ta, Ta, 9a + 2nn) = I . (IV.41) 

neZ 

Here we use the definitions 

01 (tt) = 5(rbCos9b — raCos9a — u^) 5(rbsm9b — raSm9a — u^) (IV. 42) 
for u — (u^, u^) in M^, and 

/= / Vz-e-''^°^''^/'(l)i{4ia))- (IV.43) 

This integral is nothing else than a point-to-point transition amplitude in 
cartesian coordinates. We give now a derivation of the well-known formula 
for this amplitude, by relying on our methods. 

The integrand in equation (IV.43) depends only on z^(ta), z'^{ta), hence 
as explained in paragraph A. 2. 2, we introduce the linear map L : — > 
mapping z into z{ta)- The image of the integrator T>z ■ e~^^°^^^/^ on is 
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a Gaussian integrator on R^; we need only the value of the corresponding 
covariance matrix 

^«^= [ Vz-e-'^^'^^^^/' z'^{ta)z^{ta). (IV.44) 



According to formulas (A. 57) and (A. 62), we obtain 

^aP ^ ^G-f (ta,ta) = ^S-^ih - to) . (IV.45) 



According to formulas (A. 27) and (A. 32), we transform the integral (IV. 43) 
into 





1 |2 


TT 





The function 0i [u) is a 5- factor and the final result is 
I = {s [th - ta)) exp — cxp ' 



= {s {tb - ta)) ^ exp [ ^ 



^ t}} ta J V ^ 

2\ ' (IV.47) 



S tb-ta 



iFrom the equations (IV. 41) and (IV.47) we cannot derive directly the 
point-to-point transition amplitudes in polar coordinates. We defer to para- 
graph IV. 7 a further discussion of this point. 

3. Frame bundles over Riemannian manifolds.^^ 

We consider a Riemannian (or pseudo-Riemannian) manifold M of di- 
mension d, with metric g, and the orthonormal frame bundle = 0{M) 
over M, with projection tt : 0{M) — > M. We want to choose vector fields 
■ ■ ■ , X(rf) on 0(M) such that the equation on the bundle 0{M) 



— vE'(tb,P6) = ^/i"^X(,)X(^)Vl/(t5,pfa) , (IV.48) 



For a discussion of functional integrals when paths take their values in a 
Riemannian space, see [5], [14] and [18]. 



39 



where pf, belongs to 0{M) and the first-order difi^erential operators are given 
by X(a) = ^{a)(Pb)'^i gives by projection an equation on the base space M, 
namely: 

d s 

—iP{tb,Xb) = —g^^DxD^iP{tb,Xb) 

dtb 47r (IV.49) 

= -^^^{tb,Xb) . 

Here Dx is the covariant derivative defined by the Riemannian connection, 
and "i/ — ip o 7t; moreover A is the Laplace-Beltrami operator on M. 

It has been shown in [13] that the covariant Laplacian A at a point Xb 
of M can be lifted to a sum of products of Lie derivatives h'^^X(^a)X{p) at 
the frame pb in 0{M); the integral curves of the set of vector fields |X(q,)} 
starting from pb at time tb are the horizontal lifts of a set of geodesies at 
Xb, tangent to the basis {cq,} of T^^M corresponding to the frame pb- The 
constant matrix {hap) has been chosen with the same signature as the metric 
g on M, and ^(cq, ep) = haf3- 

An explicit construction of (X^q,)} goes as follows. Let p{t) be the hori- 
zontal lift of a path x{t) in M defined by the Riemannian connection map"*^^ 
a : 0{M) L{TM,TO{M)); it satisfies the differential equation 

p{t) ^ a{p{t)) ■ xit) , p{tb)^Pb. (IV.50) 



If we put p{t) = {x{t), u{t)) and pb = {xb, Ub), a solution p{t) of the previous 
equation corresponds to the frame u{t) obtained by parallel transport of Ub, 
along the path x from Xb to x{t). The frame u{t) is also an admissible map 

u{t) : M'^ ^ T^(t)M 

i.e. u{t) maps a d-tuple into a vector whose components in u{t) are the 
chosen d-tuple; equivalently (u{t)~^x{t)) is the a-coordinate of x{t) in the 
u{t) frame. Set 

z{t) := u{t)-^x{t) = z''{t)ea , (IV.51) 

where {cq} denotes the canonical basis of the model space M^. Then we can 
express (IV.50) in the canonical form (11.21) 

p(i) = X(«)(p(t))i"(t) , p{tb) = Pb (IV.52) 

where -^(a) is defined by 

X{a){p) = cr{p) ■ ea = (o-(p) ou)-ea. (IV.53) 

More explicitly, for a given frame p at a point x of M, cr(p) is a linear 
map from T^M to TpO{M), mapping u into cr(p) ■ u ior u in T^M . 
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Here a; is a point of M and p = {x, u) a frame, where tt : — > T^M is an 
admissible map. 

If z{t) = e-oti i.e. z^{t) = 5^, the coordinates of x{t) are constant in 
the frame u{t) parallel transported along x{t). Therefore x{t) is the geodesic 
defined by 

x{tb) = Xb , x{tt) = Ca . (IV.54) 
With z(t) = Ca, equation (IV. 52) reads 

P{a){t) = X(«)(p(/3)(t)) 5f = X^a){P{a){t)) ■ 

The horizontal lift P[a.)if) of the geodesic (IV.54) is, as desired, the integral 
curve of X(^^) going through pb at time tb- With ^(a) defined by (IV. 53), p 
can be expressed in terms of the Cartan development map 

(tt o p) {t) = x{t) = (Dev z) {t) . (IV.55) 

The Cartan development is a bijection from a space of pointed paths (paths 
with a fixed end point) on Tx^M (identified to via the frame pb) into a 
space of pointed paths on M - or vice versa. Here 

Dev : VoT:c,M ^ Vx,M by z^x. (IV.56) 

The path x is said to be the development of z, if x{t) parallel transported 
along X from x{t) to Xb is equal to z{t) trivially transported to the origin of 
Ta;,M, for every t e T. 

The path integral solution of (IV. 48) is 

^{tb,Pb)^ [ Vz-exp(--Qo{z))^p{ta,z)) (IV.57) 

with p(t, z) solution of (IV.52). The path integral solution of (IV.49) is 

^{tb,Xb)= f Vz-exp(--Qo{z))(j){{Beyz){ta)). (IV.58) 

If a scalar potential is desired in (IV.48) or (IV.49), one can proceed along 
either of the methods outlined in paragraph II. 4. 

The semiclassical approximation [5] of with or without scalar po- 
tential, is considerably more complicated to compute than the semiclassical 
approximation of ^ given in section III. For an initial wave function of type 
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(III. 2) the blue-print given in Section III. 2 is complete. If one wishes to com- 
pute the point-to-point propagator on a Riemannian space M, one chooses 
the initial wave function on M to be 

<j>{x) = 5,S^) . (IV.59) 

The detailed calculation can be found in reference [5, pp. 309-311]. The devel- 
opment map cannot parametrize spaces of paths with two fixed points-*^^ but 
(IV.59) with {T)ev z){ta) substituted to x restricts the domain of integration 
appropriately. 

The following remarks simplify the calculations. 

Remark 1. According to the formulas 

I dtha,3z''{t)zf^{t)^ I dtg^,{x{t))x^{t)x''{t) 
•^"^ (IV.60) 

dtg^,{x{t))p^{t)p''{t), 

development map and horizontal lift preserve lengths and angles. 

Remark 2. If z develops into a classical paths Xd, the determinant of the 
derivative mapping Dev {z) is unity. For the proof see, for instance, reference 
[5], p. 308. 

4. A multiply connected manifold. 

The domain of integration, a space of pointed paths VxN, is the union 
of disjoint sets made of paths in different homotopy classes. We recall in 
paragraphs a) and b) earlier calculations of propagators on multiply connected 
spaces. Then, in paragraph c), we explain how these methods fit into our 
general framework. 

a) It was shown [3] in 1971 that, for a system with a multiply connected 
configuration space N, the propagator is a linear combination of propaga- 
tors 



\K 



J2 X''{9a)K^a) 
5c.e7ri(Ar) 



(IV.61) 



Each -f^(Q;) is obtained by summing over paths in the same homotopy class, say 
a; the set {x'^iQa} forms a representation, labeled A, of the fundamental 



16 Pqj. instance, two geodesies on S'^ intersect at two antipodal points; they 
are the developments of two halflines with one common origin. 
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group 7ri(A^). Since a homotopy class cannot be identified uniquely with an 
element of 7ri(A^), the propagator is defined modulo an overall phase factor. 
There are as many propagators K"^ as there are inequivalent representations 
of 7ri(iV). 

The proof of (IV. 61) uses two facts: 

i) the superposition principle of quantum propagators implies the linear 
combination of partial propagators; 

ii) the fundamental group based at a point is isomorphic to the funda- 
mental group based at another point, but not canonically so. Therefore the 
pairing ((/q,, is done by choosing an homotopy mesh (choosing a point 
for the fundamental group, and pairing one group element with one homotopy 
class), then requiring that the result be independent of the homotopy mesh. 

b) Later on [14, see also 18, p. 65] the same result (worked out for a 
different example) was obtained from stochastic processes on fibre bundles. 
The basic steps are as foUows^^. 

i) A universal covering TV is a principal G-bundle with projection 
n : N ^ N, where N = N/G and G is a discrete group of automorphisms of 
N isomorphic to the fundamental group of N. For example, N = = M/Z 
and iV = IR is a Z-principal bundle over S^. 

ii) The wave function for a system with configuration space AT is a sec- 
tion of a bundle weaMy associated to A^. i.e. a bundle whose typical fibre is 
associated to a not necessarily faithful representation of G, hence to a repre- 
sentation of a group homomorphic to G, says Go- For example ii N — S^, 
a vector bundle over with structure group ^7(1) is weakly associated to a 
Z-principal bundle over by a homomorphism of Z into U{1) mapping n 
into 

iii) There is a unique connection on N: the horizontal lift x of a pointed 
path X with fixed point x{to) = xq is uniquely defined by the lift xq of a^o- 
This unique connection defines the parallel transport of the wave function 
4){x{t)) back to a^o- 

Consider for instance the case where N — M., N = and Il{x) = e*^. 
Take for a section of a t/(l)-bundle defined by the homomorphism ha as 
above. Let x{t) be a map into S^, lifted to a map xo{t) into IR in such a way 

that x{t) = e*^°^*). The other liftings are given by 

Xk{t) = xo{t) + 27rk (IV.62) 



We refer the reader to paragraph IV. 7 for an explicit example where we 
use this strategy. 
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for k in Z. The parallel transport of (j){xk{t)) to xq is given by the formulas 

ryi^t)) = K{xo)K{xk{t)r^ Hx{t)) (IV.63) 

and 

haixkiW = haiMt))'' h^{k)-^ = e-'^°^K{xo{t)) . (IV.64) 

The map hoi is the bundle map from the Z-bundle M over to the U{1)- 
bundle corresponding to the map /Iq, : Z ^(1)- 

iv) Summing the wave function Tf^(f){xk{t)) over all /c's gives the solution 
Xq) of a parabolic equation with initial value ^{to, xq) — 4>{xo). 

In both computations outlined above, one chooses a representation of the 
fundamental group, {x^{9a)}j or the homomorphism h : G ^ Gq. 

c) We specialize our basic formula (II. 1) 

^{h, Xb) ^ [ Vz- e-^Q°^'^/'<i){xb ■ E{ta, z)) (IV.65) 

to the case where N is multiply-connected. To calculate the point-to-point 
transition amplitudes, we select (f){x) of the form S^^ (x) with a 5-factor cen- 
tered at a point Xa of N. Denote the evaluation map taking z into Xb • z) 
by £ : Zft — > A^. Since Z5 is contractible, we can lift £ to a map e : ^ N 
into the universal covering N of N and hence £ = 11 o £ . In the path integral 
(IV.65), the domain of integration is restricted by the (5-factor S^^ix) to the 
inverse image e~^{xa)- It consists of paths such that Xf, ■ T,{ta, z) = Xa and 
splits as the union of domains Z-p = e~^{xr) where xr runs over the various 
points of N mapping to Xa by 11. The labels F correspond to the various ho- 
motopy classes of paths x : T ^ N such that x{ta) = Xa, xith) = Xb- Hence 
the integral (IV.65) splits into a sum of integrals over the various subdomains 

*(tfa, Xb) = Y, [ Vz- e"-'3o(.)/s5^_^ . z)) . (IV.66) 

This equation is the justification of the heuristic idea used in (IV. 61) 
that the building blocks of K are the propagators i^(a) obtained by summing 
over paths in the same homotopy class. To explain the coefficients x'^ida) 
we can proceed as follows: using the previous notations N ^ 11, G, consider 
a homomorphism of G into U{1). The corresponding wave functions are 
functions (f) on N such that (j){xg) — X^i9)4'{x) for ^ in G and x in N. We 
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denote by z) the lifting of z) to N and generalize equation (IV. 65) 

by ^ 

^{h, Xb) = [ Vz- e-'^^°(^)/>(x6 • z)) . (IV.67) 

The function \1/ will satisfy the same transformation property as (p and when 
is a 5-factor, we can split the integration domain into subdomains Zr as 
above. We recover then the results derived in paragraphs IV 4a and b. 

5. Gauge fields. 

We begin with the case of an abelian gauge group. In physical terms, 
we consider a particle of mass m and electric charge e moving under the 
influence of a magnetic potential A, with components Aa{x) at the point x. 
We consider generally a d-dimensional space in cartesian coordinates x'^ 
(a e {1, • • • , d}) and metric |a;|^ = 5af3X°'x^. The classical Lagrangian is given 

by 

Tjfx. x) = 

2 



L{x,x) = ^\x\'^ + eAa,{x)x" , (IV.68) 



hence the action functional 

|2 



S{x) = [ dt L{x,x) = ^ / + e / Aadx"" . (IV.69) 



The equation of motion can be derived from this Lagrangian, and can be put 
into the Hamiltonian form with the following definitions 

Pa = mi" + eAa , H=\p- eAf /2m . (IV.70) 

The corresponding Schrodinger equation is obtained in the standard way by 
replacing by the operator j in the definition of H, and reads as 

^^^4-^y.(-^- -^o.{x)\ ' ^ . (IV.71) 

Our goal in this paragraph is to fit the well-known path integral solution of 
this equation into our general framework. 

It has long been recognized [19] that it is desirable to treat as a section 
of a complex line bundle (here over M*^) associated to a principal [/(l)-bundle 
via the canonical representation of acting on C by multiplication. We 
describe the main steps of this construction. 
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5.1. The invariant formalism. 

The base space is M = W^. The gauge group G is the set U{1) of complex 
numbers of modulus one g = e*®. We consider a principal bundle P, with G 
acting from the right via (p, g) ^ p ■ g, and projection 11 : P — > M. 

The connection is a differential form a; on P with the transformation 

rule 

io{p-g)=u;{p)+g-'dg. (IV.72) 

With the angular coordinate such that g = e*®, one obtains g~^dg = idQ 
for the invariant differential form on U{1), hence the Lie algebra g of U{1) 
is naturally identified with the set of pure imaginary numbers, and u is pure 
imaginary. 

For any path x : T M and any point pb in P with Il{pb) = x{th), the 
horizontal lifting is a curve ^ : T — > P satisfying the following conditions: 

mt) = xit) , atb)=Pb (IV.73) 

{u;^it),m)=0, (IV.74) 

where a^^(t) is the value of u> at the point ^(t) of P, that is a linear form on 
the tangent space T^(t)P. 

Let L be the associated line bundle. For any point x in M, a point p of 
P with n(p) = a; corresponds to an admissible map p '. C — ^ Lj; where Lx is 
the fiber of L above the point x. If is a section of L over M, its value at 
point X is an element of L^, hence p~^{il){x)) is a complex number 
In this way (see e.g. [20], vol. I, p. 404), we identify the section -0 of L to a 
function \E' : P — > C with the transformation rule 

^ip-g) = g-^-np) (IV.75) 

for p in P and g in G. 

5.2. Fixing the gauge. 

A fixing of the gauge corresponds to a section s : M ^ P of the principal 
bundle. When s is chosen, we may identify P with M x G in such a way that 
p = s{x) ■ g in P correspond to {x, g) in M xG. A section ijj of the line bundle 
L corresponds now to a wave function, that is to a complex- valued function on 
M and the corresponding function on M x G is given by ^{x, g) = g~^'ip{x) 
and conversely i(^{x) = ^{x,l), or intrinsically i/j = ^ o s. 

The differential form cj on P gives by pull-back via s : M — > P a purely 
imaginary differential form on M, to be written as —^A to fit with standard 
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physical dimensions. Hence the differential form ^4 on M can be written 
as Aci{x)dx°' and the functions Aq,(x) are the components of the magnetic 
potential. On M x G the differential form uj is given by 

u = g-^dg -^-^A^i(de- ^Ac^dx"^ (IV.76) 

(for g = e*®). 

Let X : T — s> M be a path. The horizontal lifting ^ of a; is now described 
by ^{t) = (ic(t), e*®*-*-*) and since uj induces a zero form on the image ^(T) C 
M X G, we obtain the differential equation 



e=^A^{xit))x-{t) (IV.77) 



(see e.g. [18, pp. 64-65]). 



Changing the gauge corresponds to choosing another section si : M — > P. 
There exists then a function i? : M — > M such that si{x) = s{x) ■ e~*^^(^)/^. 
In the new gauge, the section of the line bundle L corresponds to a new wave 
function 

iPiix) = e'''^^'^^/^^{x) (IV.78) 
and the new components of the magnetic potential are given by 

Al{x)=A^{x) + ^-^. (IV.79) 



5.3. Path integrals. 

We revert to the notations in paragraph IV. 1. The constant A is again 
{h/mYf^ and we parametrize the paths x with x{t},) = Xb by 

x{t, z)=Xb + \z{t) (IV.80) 

where z runs over the space Z;,. The horizontal lift of the previous path is 
given by 

e(t,^) = (x(t,^),e'«(*'^)) . (IV.81) 

Taking into account the differential equation (IV.77), we obtain the following 
differential system 

dx'^it) = Xdz'^it) for a in {I,-- -,4 

(IV.82) 

dQ(t) = fA^{x{t))dz^{t). 
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This system has the canonical form (11.21) where the vector fields ^(i), • • • , -^(d); 
Y are given by 

Notice that for ^{x,Q) = e~'^®ip{x), we have 

£x(„)^'(x, 6) = e-*® ■ XDai^{x) (IV.84) 
with the differential operator (see e.g. [20, p. 405]): 

Da = ^-jAM- (IV.85) 



Our general partial differential equation 



^ = f * (IV.86) 



a. 



translates now as 

a 

for the ^/'-compoment of For s = i, this equation coincides with the 
Schrodinger equation (IV. 71). 

Feynman path integral solution to this equation reads as follows: 

i^{tb,Xb)= [ Pa; • e^^(^)/V(x(ta)) (IV.88) 

where the action S{x) is given by equation (IV.69). Parametrizing the paths x 
in Vb by the paths z in Zb (see equation (IV. 80)), we can rewrite the previous 
path integral as^^ 

il^{tb,Xb) = Vz-e^'^^^^'^expi^j j^A^{x{t,z))dx''{t,z)^ (l>{x{ta, z)) . 

(IV.89) 



18 We give the formula in the oscillatory case s — i. The reader is invited 
to work out the formulas for the case s — 1. The definition of Qo{z) is given 
in equation (IV.22). 
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Replacing ^l>{x) by ^{x,Q) = e and similarly 0(a;) by ^{x,Q) = 

e~*®0(a;), we can absorb the phase factor and obtain 

Xb, Qb) = f Vz- e"^^°(^)$((a;6, Qb) ■ S(ta, z)) . (IV.90) 

The transformation 2;) of the bundle space M x G takes {xb,Qb) 
into {xb + Xz{t), Qb — ^ J-i-'' Act {xb + Xz{t)) dz"{t)) and corresponds to the 
integration of the differential system (IV. 82). 

5.4. Various generalizations. 

a) It is easy to incorporate an electric potential V. The complete action 
functional is now 

S{x) = ^ / + e f A^dx'^ - Vdt . (IV.91) 

Feynman solution (IV. 88) is still valid and can be made explicit as 

i^{tb,xb) = Vz-e'^^'^^^'^exi^l^j j^A^dx'^-Vdt^4>{x{ta,z)) (IV.92) 

where the line integral Jrj.Aadx'^ — Vdt is calculated along the path x{-,z). 
The Schrodinger equation reads as follows: 



b) A non-abelian gauge group is for instance G = U{N); a more general 
compact gauge group can always be realized as a closed subgroup of some 
unitary group U{N). We mention a few of the required changes. The Lie 
algebra g is the set of matrices of the form lE where E is an N-hy-N hermitian 
matrix. The connection form u on the principal bundle P takes its values in 
Q and the transformation rule is now 

w(p ■ 9) = 9~^^{p)g + g'^dg (IV.94) 

where 00 {p), g, dg are N-hy-N matrices and the products are given by matrix 
multiplication. 

For the associated vector bundle L, we consider the natural action of 
U{N) on the complex vector space C^. Hence a section ip oi L corresponds 
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to a function : P — > such that '^{p ■ g) = g ^ ■ '^{p) for p'm. P and g in 
G=U{N). 

In a given gauge the form uj is given by uj{x, g) = g~^dg — ^A{x)^ where 
A{x) = Aa{x)dx°' is a hermitian N-hj-N matrix of differential forms on M. 
A gauge transformation is given by the formulas 

iji{x) = U{x)-^ij{x) (IV.95) 

Al{x) = U{x)-'Aa^{x)U{x) + '-^U{x)-'^U{x), (IV.96) 

e ox"' 

where U{x) is a unitary matrix depending on the point x of M. 

The horizontal lift of a path a; : T ^ M is of the form ^{t) = {x{t), U{t)) 
where the unitary matrix U (t) satisfies the differential equation 

U{t) = -U{t) ■ Aa{x{t))x''{t) . (IV.97) 

We can solve this equation in the standard way using time-ordered exponen- 
tials Texp. 

The Schrodinger equation is still written in the form (IV. 87), but Dq, 
is now a matrix of differential operators, namely ■ 1 — ^Aa{x) where 
Aa{x) is an N-hy-N hermitian matrix. In the path integral (IV.89) replace 
the exponential factor by 

Te^p(^^ J^" Aaixb + Xz{t))dz''{t)^ . (IV.98) 



c) We could consider gauge groups over a curved manifold and combine 
the results of paragraph IV. 3 with those of the present paragraph. 

6. A symplectic manifold. 

Let N he a symplectic manifold of dimension d = 2n. The mani- 
fold A4 represents the classical phase space of a physical system - usually 
the cotangent bundle T*Q of a configuration space Q but not necessarily. 
Paths in M. have n initial and n final boundary conditions. These boundary 
conditions are consistent with the requirement of quantum uncertainty, and 
they imply a choice of polarization^^. Hence, it is possible to cover M. with a 

Contributed by John LaChapelle. 

Roughly speaking, a polarization is a foliation of M. whose leaves are 
Lagrangian submanifolds of dimension n. 
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family of open subsets {Ui} such that each Ui is diffeomorphic to a product 
of two transverse Lagrangian submanifolds Li x L'-. For simplicity consider a 
symplectic manifold which admits global transverse Lagrangian submanifolds 
L and L', and identify M. with L x L'. Set x{t) in the form (^^(t), a;L'(t)) 
with XL{t) in L and a:L'(^) in L', in such a way that XL{tb) = and that 

XL'{ta) = Xa- 

We require the group generated by the transformations S(T, z) to be a 
subgroup of the group of symplectomorphisms which leave the polarization 
invariant. This implies that the set is of the form \^X(^a)i ^(a')} such 

that X(^a){xL'{t)) = 0, X(„/)(xL(t)) = 0, and = 0. Here a E 

{l,---,k}, a' G {/c + 1, ■ ■ ■ , d}, and /c is a fixed integer between 1 and d, 
possibly, but not necessarily, equal to d/2. Consequently, a path satisfies the 
differential equations 



dXL{t) = X^a){xL{t)) dz^ + Y{xL{t)) dt 
dXL' (t) = {xL' (t)) dz''' + Y{XL' {t)) dt , 

and the general formula (II. 1) becomes 

X(l){xb ■ T,{ta, Zl) , Xa ■ 11(^6, Zl')) , 



(IV.99) 



(IV.lOO) 



where now /i"^ ~ ( h'^'^' ) ' ^^^^ space of paths : T ^ Ml^ 

such that z^itiy) = 0, and Z/,/ is the space of paths z^i : T such that 

ZL'ita) = 0. Each functional integral separately satisfies a partial differential 
equation: 

(IV.lOl) 



where *l := *|r and '^l' := *| 



L" 
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6.1. The case of a cotangent bundle. 

Choosing an initial function (p is choosing a transition amphtude. We 
consider the case where Ai is the cotangent bundle of a flat configuration space 
Q, hence M. — Q x P where Q and P are finite-dimensional vector spaces 
in duality. In order to define position-to-position transition amplitudes., we 
choose the initial function of the form (l){q,p) = 5{q — Qa). Equation (IV. 100) 
yields 

}C{qb,tb;qa,ta) = / VzqVzp ■ ex.p(--Qo{zQ, zp)) h{zQ, zp) 
JzqJzp ^ s ) 

(IV.102) 

where the integrand is given by 

h{zQ, Zp) ^5{qb- zq) - qa) . (IV.103) 

In order for the transition amplitudes to be consistent with the initial 
wave function, we require lim )C{qb, tb', qa, to) to be equal to 5{qb — qa)- But 

the integrand h{zQ, zp) tends to 6{qb ■ T,{tb, zq) - qa) = 5{qb - qa) when ta 
tends to tb, a limit independent of zg, zp. Hence, provided we can interchange 
limit and integration, we get 

lim K:{qb,tb;qa,ta) = / VzqVzp ■ ex.p(--Qo{zQ, Zp)) 5{qb - qa) ■ 

ta^h JZp ^ S / 

Hence this is equal to S{qb — qa) by the normalization of our integrator. 
We handle the other cases in a similar way. 

a) Momentum-to-position amplitude ]C{qb, tb'iPai ta)' use the initial func- 
tion /i-"'/2e*«i''P/'^ and the integrand 

h-^/^exp(^pb ■ {Pa ■ S(tb,Zp))^ . (IV.104a) 

b) Position-to-momentum amplitude lC{pa, ta', qb, tb)'- use the initial func- 
tion /i-'^/2e~*9 Pa./'i and the integrand 

/i-"/2exp(^-^ {qb ■ i:{ta,ZQ))-pa^ . (IV.104b) 

c) Momentum-to-momentum amplitude )C{pb,tb',Pa,ta)' use the initial 
function 5{p — pb) and the integrand 

5{pa-T.{tb,Zp)-pb) . (IV.104e) 
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It follows easily from these definitions that 

K:{pb,tb;Pa,ta) = h-^/^ [ dqk e-''^''-P''/''IC{qb,tb;Pa,ta) (IV.105) 

Jo 



'Q 

K:{qb, tb-, Qa, to) = /i^"/' ^ dpb e''"'-P'/''K:{pb, tb-, qa, ta) . (IV.106) 

6.2. Coherent states. 

More general transition amplitudes, which cannot be interpreted as posi- 
tion-to-momentum transitions, are possible by choosing more complicated 
initial wave functions, different polarizations, and/or by having non-trivial 
phase spaces. For instance, coherent state transitions can be calculated - 
given a suitable characterization of coherent states. 

Choose a Kahler polarization on a non-trivial symplectic manifold, so 
that is a Kahler manifold M of complex dimension n. It is convenient 
(though not necessary) to take the phase space to be the product manifold 
M X M with coordinates (C;^) and complex dimension 2n. Following the 
work of Berezin [21] and Bar-Moshe and Marinov [22], we use generalized 
coherent state wave functions 0^'(C) — expi^(C, c') where K{(,(') is the 
Kahler potential. Consider the function : M x M ^ C given by C ) = 
g4c.Ca) and take : T ^ C" such that zmitb) = and zj^ : T ^ such 
that z-j^{ta) = 0. Then 

JZm -'z— 

are generalized coherent state transition amplitudes. Note that }C{Cbj tb'tCaj to) 



= /C(C6,4;Ca.^a), and lini /C(C6, tb; Ca, ^a) = expK{CbXa) provided the 
limit can be taken inside the integral. 

7. A simple model of the Bohm-Aharonov effect. 

We illustrate on an elementary example various techniques expounded in 
this section. Our basic manifold is the euclidean plane with the origin 
removed. As in paragraph IV.2, we denote the cartesian coordinates by x^, x^. 
The universal covering of N is the^set N of pairs of real numbers r, 9 with 
r > 0, and the covering map 11 : iV — > A?" is described by equation (IV. 31), 
namely = rcosO, x^ = rsinO. 

We introduce a magnetic potential A with components 

A, = —^ , A2=-^. (IV.107) 

2'k\x\ 2'k\x\ 
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The magnetic field has one component B perpendicular to our plane, given 

and an obvious calculation gives B = outside of the origin. But the circu- 
lation ^ Aidx^ + A2dx^ on any loop around the origin is equal to F, hence 

B^F5{x). (IV. 109) 

Physically, we have a wire perpendicular to the plane through its origin, car- 
rying a magnetic flux equal to F. 

Denote by m the mass and by e the electric charge of a particle. For a 
path X : T ^ N, with components x^{t), x^{t) at time t, the action is given 

by 

S{x) ^ So{x) + Sm{x) , (IV.llO) 

with the kinetic action 
and the magnetic action 

We fix a point Xh in A?", and denote by X5 the set of paths a; : T — > AT for which 
the kinetic action Sq{x) is finite (that is the L^'^ paths) with the boundary 
condition x{tb) = 2^5. On X;,, we denote by Vx the translation-invariant 
integrator normalized by 

j Vx ■ ey.^{^So{x)^ ^ I . (IV.113) 
The free position-to-position transition amplitudes are defined by 

{tb,Xb\ta,Xa)Q= J Vx ■ exp(^^So{x)^ S{x{ta) - Xa) . (IV. 114) 

According to the calculation made at the end of paragraph IV. 2, we get^^ 

— TTii I i \ 

{tb,Xb I ta,Xa)Q = ^^^^ _ ^ ^ exp\-So{xc\) j , (IV.115) 



We express our formulas directly in terms of the path x, and have no 
need for the scaling x = Xz introduced earlier. 
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where Xd is the classical path of a free particle, namely: 



tb - ta 

Hence, we obtain explicitely 

Z tb — ta 

that is, the WKB approximation (IV. 30) is exact in our case. 

We consider now the transition amplitudes in the given magnetic poten- 
tial: 

(tb, Xb I ta, Xa)p = ^ • {j^S{x)^ S{x{ta) - Xa) ■ (IV.118) 

According to (IV. 87), we consider the Schrodinger equation 



dt 2m 

d 



{Dl + Dl)ijF (IV.119) 



with the differential operators = — ^^a- Explicitly, we obtain 



dt 2m 

with 



LV'F, (IV. 120) 



d \^ f d \^ 2ci f I d 



The constant c is equal to eF/h; it is dimensionless. The solution to the 
equation (IV.119) is given by 

1pF{tb,Xb)^ dXa{tb,Xb\ta,Xa)p'lpF{ta,Xa) , (IV.122) 

Jn 

a formula equivalent to our more familiar one 

i^Fitb,Xb) = ^ • e^p(^^S{x)^ (l>{x{ta)) (IV.123) 

if we take into account the initial wave function (p{x) — il^Fita, x) at time ta- 

55 



Wc lift now everything to the universal covering A^.^ By means of the 
formulas (IV.31), the wave function is now a function iljF{t,r,9) with the 
restriction 

ijFit,r,6 + 2n)=i^F{t,r,6) (IV.124) 

that is, a function on T x invariant under the group Z acting on by 
{r,9) ■ n = {r,9 + 27rn). The Schrodinger equation (IV. 120) keeps its form 
with L changed into the new operator 

~ 1 Id 2ci d c2 

This operator can be written as L = e*^*^ Ae~'^*^ where A is the Laplacian in 
polar coordinates, namely: 

1 Id 

Hence the new wave function 

Mt, r, 9) = e-'^'^ipFit, r, 9) (IV.127) 

satisfies the equation 

that is the Schrodinger equation for a free particle, written in polar coor- 
dinates. As we shall see in a moment, equation (IV.127) expresses a gauge 
transformation which apparently removes the magnetic potential. But the 
periodicity condition (IV.124) transforms into 

V^o(i, r,9 + 27r) = e-^'^'^Mt, r, 9) , (IV. 129) 

hence the wave function •0o on N is not the lifting of a function •0o on N . 
Rather, the linear representation n i— > e^'^*"'^ of the group Z into U{1) defines 
a line bundle on N = ]R^\{0} and i/jq corresponds to a section of this 
line bundle. Otherwise stated, despite the fact that the magnetic field B = 
dAi/dx^—dA2/dx^ is identically zero on N , there is no single- valued function 
R on N such that A^y — —dR/dx"": the space N is multiply connected. 

To explain gauge transformations, let us introduce the principal bundles 
P = N X U{1) and P = N x U{1). According to the general formulas the 
connection form lvf on P is given by 

ajF = i(dQ-^A) =i(de- c ^'"^^' ~ ^'"^f | , (IV.130) 
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where is the angular coordinate on U{1) (taken modulo 2n). When there 
is no magnetic field, it reduces to coq = idO. We cannot transform ti;i7;^into 
o^o by a gauge transformation, but if we lift these differential forms to N, we 
obtain 

u;f = iide - cdO) , ujq = idQ . (IV.131) 

The transformation U oi P into P taking (r, 9, O) into (r, ^, G + c9) is an 
automorphism of ?7(l)-bundle, and U* but U is not the lifting of an 

automorphism U of P, except when c is an integer. 

To the wave function '4>f{x) on N we associate the function 

^r^(x, 6) = e-'^il^F^x) (IV.132) 

on P. The Schrodinger equation (IV. 119) translates into our standard form 

^ = ^h'^PCxP *F (IV.133) 

provided we take /Iq,/? equal to mSafs/h, with the vector fields 

d cx^ d 



'(1) dx" bl'^ 



'1) ux \x\ 

(IV.134) 

r — o j_ cx o 

We lift now everything to P, so ipp lifts to V'f, *f to and the lifted vector 
fields are given by^^ 



_ /3 . 9 sin 6* . ( d ^ d \ 

- cosy -g^-^r- ym^^'m) 

(IV. 135) 

L -sm9 ^ + ^.^^ + c-^j . 



(2) 

The Schrodinger equation (IV.133) remains the same with replaced by 
*F and X[^^ hyXf^y 



We describe the effect of the gauge transformation U given by U{r, 9, 0) = 
(r, 6*, + c9). We noticed that U*uf is equal to ujq, and similarly, we get 
"ifFoU = \E'o where \E'o is defined by analogy to namely 



*o(r, 9, 0) = e-'®Mr, 0) . (IV.136) 



22 Yor F — 0, they reduce to the vector fields given in equations (IV. 36) 
and (IV.37). 
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Dually, U transforms the vector field into X^^^ for a = 1 or 2. Hence 

is a solution of the Schrodinger equation (IV. 133) if, and only if, is a 
solution of the corresponding equation for F = 0. 

We conclude this paragraph by a discussion of path integrals. We parame- 
trize paths in X;, by paths in Z5, using the correspondence 

x(t) = Xb + z(t) . (IV.137) 

Fix and 6*6 in such a way that 

Xb= {rb cos 9b, Tb sin 9b) . (IV. 138) 

Then x can be lifted in a unique way to a path x : T ^ N of the form 

x{t) — (r(t), 9{t)) such that x{tb) = (r;,, 9b), that is r{tb) = Vb and 9{tb) = 9b. 
The initial point Xa = x{ta) is lifted to x{ta) = {ra,9a) where r{ta) — r^, 
9{ta) = 9a- Moreover, the magnetic action of the original path is given by 

SM{x)/h = c-{9b-9a) . (IV.139) 

Consider the space Xo,6 of paths x :T ^ N such that x{ta) — Xa, x{tb) — Xb- 
Then two^such paths are in the same homotopy class if, and only if, the lifted 
paths in N correspond to the same determination of the angular coordinate 9a 
of Xa- An equivalent condition is that they have the same magnetic action- 
Using the connection of the principal bundle P over A^, we can define the 
horizontal lifting of x; it is the unique path of the form y{t) = (r(t), 9{t), ^{t)) 
on which induces a zero form, that is Q{t) — c9{t) is a constant in time. 
For the phase factors, we get 

^i&{tb) — ^'ie{ta)^iSM{x)/h ^ (IV. 140) 

Moreover, the horizontal lifting is a solution of our standard differential equa- 
tion dy = X(a)(y) • dz"^. 

^From all this, it follows that the transition amplitudes as defined in the 
Feynman way (IV. 118) agree with our standard expression (see for instance 
(IV.38)). The result obtained in paragraph IV. 2 (see formulas (IV. 41) and 
(IV. 47)) can be restated as 

{tb, Xb I ta, Xa)^ = — XI (*^' I ^« + ^"'^) • (IV.141) 

A similar reasoning yields a more general formula: 

{tb, Xb I ta, Xa)^ = -Yl e^^^^^-^'^-^"-) {tb, rb, 9b I ta, ra, 9a + 2n7r> . 



(IV. 142) 
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We can invert this formula^^ and obtain 

{h,n,9b I ta^Ta^Oa) = r„ ^ e-^'^^^""^") (tft, Xft | ta^Xa)^^^^ • (IV.143) 

Hence, from the knowledge of the transition amplitude (tb^xi, | ta,Xa)p as 
a function of the magnetic flux F , we can infer the value of the transition 
amplitudes in polar coordinates. 



Notice that the transition probabilities are functions of F which admit a 
period h/e, a, well-known physical effect. 
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JAppendix A 
Functional integration 



In this appendix, we develop the basic properties of our integrators. In 
the oscillating case (s = i), our theory is, up to some inessential changes in 
notation, the same as the one expounded by Albeverio and H0egh-Krohn in 
[9]. The introduction of a parameter s equal to 1 or z enables us to treat in 
a unified way the oscillating integrators of Fresnel type e^^^^^^Vx and the 
Gaussian integrators of type e~'^^^^'>Vx (with Q{x) > for a; in the 
latter case). We content ourselves by giving here the basic formulas and the 
computational tools. In the Gaussian case, we would have to justify these 
formal manipulations, since, as it is well-known, the L^'^ functions are a set 
of measure for the Wiener measure, and we are considering functions on 
this set of measure 0. Our claims can be fully vindicated, but we defer the 
complete justification to another publication. 

1. Gaussian integrators. 
1.1. The setup. 

We denote by X a real, separable, Banach space, by X' its dual and by 
{x', x) (or sometimes (x', x)y^ the duality between X and X'. We suppose 
given a continuous linear map D : X ^ X' with the following properties: 

- (symmetry) {Dx, y) = {Dy, x) for x, y in X; 

- (invertibility) there exists a continuous linear map G : X' — > X inverse 
of D, that isDG=l and GD = 1. 

Out of these data one constructs two quadratic forms, Q on X and W on X', 
by the rules 

Q{x) = {Dx, x) , W{x') = {x', Gx') . (A.l) 

They are related to each other as follows: 

Q{x) = W{Dx) , W{x') = Q{Gx') (A.2) 

for a; in X and x' in X'. 

We denote by s a parameter equal to 1 or i. The function e"'^^^ on X' 
is continuous. It is bounded in the following cases: 

- s = z; 
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- s = 1 and W{x') > for 7^ in X'. Equivalently, by (A.2), s = 1 
and Q{x) > for a; 7^ in X. 

1.2. The oscillatory case {s = i). 

The integrator T>x is characterized by the following integration formula: 

/ Vx ■ exp{7riQ{x) - 27rz {x',x)) = eKp{-7riW{x')) (A.3) 

for every x' in X'. This relation should be interpreted as follows. 

Since the metric space X' is complete and separable, we know the no- 
tion of a complex bounded Borel measure^^ on X'. The Fourier- Stieltjes 
transform jFyU of n is given the customary definition: 

{J'fJ,){x)= [ d/x(a;')e-'"'(^''^). (A.4) 
ix' 

It is a continuous bounded function on X. 

We denote by^^ -^(X) the set of functions on X of the form $^ = e'^^'^ ■ 
J-'H, where runs over the measures on X'. Since the map n 1— > J^/j, is injective, 
JF(X) is a Banach space, the norm of the function $^ being taken equal to the 
total variation^^ Var(//) of /i. On this Banach space, one defines a continuous 
linear form, denoted as an integral, by 

/ Vx-^^{x)= [ dfx{x')e-'''^^'''\ (A.5) 
JX ix' 

Hence by definition, we have 

/ Vx / diJ,{x') exp(7TiQ(x) — 27ri {x' ,x)) = / diJ,{x') exp{—7riW{x')) . 
Vx Vx' Jx' 

(A.6) 



That is: a cr-additive functional from the cr-algebra of Borel subsets of 
X' into the complex numbers (see [23] for instance). We use a simplified 
terminology by refering to /i as a "measure" on X'. 

The initial stands for Fresnel or Feynman according to the worshipping 
habits of the reader. 

According to the standard definition, this is the l.u.b. of the set of num- 

p 

bers ^ |A*(^i)l where {Ai, ■ ■ ■ ,Ap) runs over the set of all partitions of X' 

i=l 

into a finite number of Borel subsets. 
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Formally, this is the equation obtained by integrating equation (A. 3) w.r.t. 
dfj,{x') and then interchanging the integrations: 



/ dn{x') Vx = Vx dn{x') . 
ix' ix ix ix' 

The space ^(X) of Feynman-integrable functionals on X is invariant 
under translations by elements of X, and so is the integral, namely: 

/ Vx ■ F{x) = I Vx-F{x + xo) (A.7) 
Jx Jx 

or, in shorthand notation, Vx = V{x + xq), for any fixed element xq of X. 
According to equations (A.l) and (A. 2), formula (A. 3) can be rewritten 

Vx ■ e'^^^^^-G^') = 1. (A.8) 



Hence assuming the invariance under translations of the integral, the normal- 
ization of Vx is achieved by 



Vx ■ e'^*^^^^ = 1. (A.9) 

X 



1.3. The positive case (s = 1). 

Henceforth, we assume that s = 1 and that Q is positive-definite (that is 
Q{x) > for X ^ in X). We take as basic integration formula 

/ Vx ■ eiip{-nQ{x) -2TTi{x',x)) ^ exp{-nW{x')) . (A.IO) 

We can interpret this relation as above. The conclusion is the integration 
formula 

/ Vx / diJ,{x')e-xp{—7TQ{x) — 27Ti{x',x))= / dfx{x') exp{—7TW{x')) . 
ix ix' Jx' 

(A.11) 

But the space .F(X) is no longer invariant under translations. To restore this 
invariance, we have to replace Fourier-Stieltjes transforms by Laplace-Stieltjes 
transforms. For that purpose, we have to consider the complex dual space 
X^ consisting of the continuous real-linear maps x' : X — > C and measures /j, 
on Xj^. We leave the details to a forthcoming publication. 
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Formulas (A. 3) and (A. 10) are obtained as the specializations of formula 
(1.7) for s = z and s = 1 respectively. 

2. Linear changes of variables. 

2.1. The case of Lebesgue integrals. 

Let Y be another real, separable, Banach space, and let L : X ^ Y be 
a Borel-measurable map. If is a measure on X, its image under L is the 
measure cuy on Y defined by (jjy{B) = (jj^(^L~^ (B)^ for any Borel subset B 
of Y. In functional terms, this definition is tantamount to 

/ dij^ix) g{L{x)) ^ f du^{y)g{y) (A.12) 

for any bounded (or non- negative) Borel-measurable function ^ on Y. 

_ Assume now that Lis linear and continuous. We consider the transpose 
L of L, that is the map L : Y' — > X' such that 

{Ly',x)x={y',Lx)^ (A.13) 

for a; in X and y' in Y'. The measures wx on X and on Y being as 
above, introduce their Fourier- St icltjes transforms JFcux and J-'uy respec- 
tively. Hence JFcux is a bounded continuous function on X', and similarly for 
J-'uy- on Y' . By definition, we have 

.Fa;x(a;') = / cZa;x(a;) exp(-27ri (x', a;)^) (A.14) 
ix 

J^io^iy') = dio^iy) exp(-27ri {y' , y)^) . (A.15) 



By specializing g{y) = ex.p{—2ni (y' , y)^) iiito formula (A.12) and taking 
into account formulas (A.13) to (A.15), we get J-'uy-{y') — J^u)y^{Ly') for every 
y' in Y'. Hence the composition formula: 

Tlo^ = Tujy^o L . (A. 16) 
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[The figure can be found in the Journal of Mathematical Physics 36 p. 2284] 



The case of a translation is similar. Assume now that xq is a given 
vector in X and denote by T the translation taking x into x + Xq in X. If u 
is any measure on X and u>xo its image under the translation T, a suitable 
specialization of formula (A. 12) gives the Fourier transform of cOxq, namely 

(JTo;^ J {x') = exp(-27rz (x', xq)) ■ Tu{x') . (A.17) 

2.2. Infinite- dimensional integrators. 

We go back to the setup of paragraph A. 1.1. We denote by Vx the 
integrator characterized by 

j Vx-eK^{^-^Q{x)^ ■e-xp{-2TTi{x',x)) = eyc£,{-TTsW{x')). (A.18) 

Formally, this means that the integrator defined by 

Vuj{x) = exp(^--Q(x)) • Vx (A.19) 
has a Fourier transform equal to exp{—7rsW), namely: 

/ Vu;{x) exp{-2ni {x' , x)) = exp {-nsW{x')) . (A.20) 

This can be interpreted as follows: the integrator Vu is a continuous linear 
form on the Banach space of Fourier-Stieltjes transforms J-'n, given by 



f Vu;{x)J^n{x)^ dii{x')exp{-7TsW{x')). (A.21) 
X ix' 



This is another form of the Parseval relation. 
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We proceed to study the image of the integrator Vuj under a hnear map. 
We use the following notations: 

- Qx is a quadratic form on X, with inverse Wx on X'; 

- Qy is a quadratic form on Y, with inverse WV on Y'; 

- L is a continuous linear map from X into Y. 

We assume that L is surjective, hence L is injective, and that the quadratic 
forms Wx and are related by 

= Wx o L . (A.22) 

Consider now the integrators^^ 

Vuj^{x) = exp(^--Qx(a^)) • T^x (A.23) 
Vu^iy) = exp(-^gY(y)) • Vy. (A.24) 

The Fourier transforms are given respectively by 

.Fa;x = exp(— TTsWx) , J^oJy = exp(— ttsWV) (A. 25) 

and according to formula (A.22), we obtain 

J^ujy = Tuiy^ o L . (A. 26) 

This is the same as formula (A. 16), hence ivy is the image of ux. under the 
linear mapping L. Explicitly stated, we obtain the integration formula: 

J^Vx ■ exp(^~Qx{x)^ ■ g{Lx) ^ J^Vy ■ expl^—Q^iy)^ ■ g{y) . (A.27) 

This holds if is a Fourier-Stieltjes transform for some measure v on Y'. 

The relationship between the quadratic forms Qx on X and Qy on Y is 
expressed by the relation (A.22) between their inverses. Let Dx : X — > X' be 
the continuous linear map such that Q~s.{x) = {D:x.x,x)^ and {D-^xi, X2)x = 
(il>x3^2, a^i)x' define similarly D^. Between the inverses Gx of Dx and 
Gy of -Dy there holds the relation: 

Gy = o Gx o Z . (A.28) 



In the space X, the integrator Vx depends on Qx and s, and should be 
written more explicitly as Vg^g^x. Similarly for Vy in space Y. 
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We mention a few particular cases: 

- if L is invertible, with inverse L~^, then Qy = Qx ° L~^; 

- if Qx is positive-definite, then Q^{y) is the infimum of Qx over the 
set of elements a; in X such that Lx = y; 

- for any given y in Y the equation 

D^x = (l o L>y) (y) (A.29) 
has a unique solution x = x{y) in X and we obtain 

Q^{y) = Q^{x{y)) . (A.30) 



3. Examples and applications. 

3.1. The finite-dimensional case. 

We record here the basic formulas. Assume that X is of finite dimension 
d; after choosing a linear frame, we represent a vector a; by a column matrix 
(x") where a G {1, ■ ■ ■ , (i}. The elements of X' are represented by row matrices 
and the duality is given by (x' , x) = x'^x" (Einstein's summation convention). 
The volume element is given hj dx — dx^ ■ ■ ■ dx"^. 

The quadratic forms Q and W correspond to symmetric matrices, namely: 

Q{x) = haox'^xf^ , W{x') = h^f^x'^x'fj (A.31) 

with hajsh^^ — 52 . The integrators are given by 

Vx=\detha0\^^^dx^---dx'^ (A.32) 
Vu}(x) = e^p{-7Thapx°'x^) ■ Vx (A.33) 

when s = 1 and the matrix {hafs) is positive-definite. In the oscillating case, 
we have to multiply Vx by e'^^^'^~P^/^ where the symmetric matrix (ha/s) has 
p positive and q negative eigenvalues. 

The quadratic form W gives the covariance. More precisely, one obtains 

/ Vu;{x) {x',xf = —W{x') (A.34) 
ix 27r 

hence 



/ Vuj{x) 
ix 



x'^x^ = ^/i"^ . (A.35) 
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3.2. Image under a linear form. 

Let x' in the dual X' of X. By specializing the results of paragraph A. 2. 2 
to the linear map L : x i— {x',x) from X into M, we obtain the following 
result. We identify M with its dual, hence L takes a number A to Xx' in X'. 
With Qx. equal to Q, hence Wx. to , we obtain 

W^{X) = X'Wix'), (A.36) 

hence 



[u) = v?lW{x'). (A.37) 
^From formula (A. 27), we deduce the following integration formula: 

/ Pa;-e-'^Q(^)/*^((a;',a;)) =C / due'''''^ I '"^^'"'^ g{u) (A.38) 

where the normalization constant C is l/(sM^(x'))^/^ (principal branch of the 
square root). More explicitly: 

- for s = 1, hence W{x') > 0, then C — ^ 



for s = z, and W{x') > 0, then C 



/4 



TTi/4 

- for s = i, and VF(a;') < 0, then C = , 

^\W{x')\ 

If we take in particular g{u) = v?, we get 

I Vx ■ e-'^Q(^)/^ {x', xf = —W{x') . (A.39) 
ix 27r 

By polarization, we obtain the more general formula 

/ Vx ■ e-'^^(^)/^ {x', x) {y', x) = — {x' , Gy') (A.40) 
ix 27r 

for x' and y' in X', thus giving the covariance of our integrator. 

Remark. In the application IV. 1 to point-to-point transitions, we need to 
restrict the domain of integration (see also paragraph A. 3. 8). This can be 
achieved by inserting a delta function 5{{x'q,x)) and then using equation 
(A.38): we restrict the domain of integration from the space X to the hyper- 
plane Xq with equation (xq, a;) = 0. 
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In general, any integral of the form 

with fixed elements x'l, ■ ■ ■ , x'^ in X' can be reduced to an n-dimensional inte- 
gral. For instance, combining formulas (A. 38) and (A. 40), the Green's func- 
tion Ga,b{t,U') on Za,b given by formulas (A. 58) and (B.28), can be evaluated 
using an integral 

J Vbz-exp(^~^j dtz{tf^5{z{ta))z{t)z{u) 

and reducing it to a 3-dimensional integral. 
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3.3. Space of paths of finite action. 

The time interval T = [ta^^b] being given, we denote by L^'^ (or more 
accurately L^'^(T)) the space of real-valued functions z{-) on T with square- 
integrable derivative i(-) and we define the quadratic form Qq on L^'^ by 

Qq{z) = J dtz{tf. (A.41) 



The definition of L^'^ can be rephrased as follows: the function z belongs to 
L^'^ if and only if there exists a function z in L^(T) such that 

z{t')-z{t) = duz{u) (A.42) 

whenever t, t' are epochs in T such that t < t'. The function z is then defined 
up to a null set; indeed by Lebesgue's derivation theorem, one gets 

z{t) =lim {z{t + r) - z{t)) /r (A.43) 

r=0 

for almost all t in T. By Cauchy-Schwarz inequality, one deduces from (A.42) 
the inequality 

\z{t') - zit)f < Qoiz) ■ \t' - t\ (A.44) 

for t, t' in T. Hence any function z in L^'^ satisfies a Lipschitz condition of 
order 1/2, and a fortiori it is a continuous function. 

The quantity Qo{z) in equation (A.41) can be calculated using a dis- 
cretization of time. Consider a subdivision T of the time interval T by epochs 

< ^0 < ^1 < • • • < ^7V-1 < ^ AT < ^6 • 



Set Ati = ti — ti-i for < i < N + 1 with the convention t_i = ta, tN+i = tb 
and denote by 5{T) the largest among the increments At^, that is the mesh 
of the subdivision T. For any function 2; : T — > R set 

Zi = z{ti) , Azi = Zi- Zi-i (A.45) 

for < i < + 1. The quadratic variation of z w.r.t. the subdivision T is 
defined as 

N 

QT{z)^Y.^Azif /Ati. (A.46) 

i=l 
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Then the function z belongs to L^'^ if and only if the set of quadratic variations 
Qt{z) is bounded when T runs over all subdivisions of T. Then 



Qo{z) = 1. u .b. Qr{z) . (A.47) 



More precisely, for any sequence of subdivisions T(n) (for n e N) whose mesh 
S{T{n)) tends to 0, one gets 



Qo{z) = \im QT(n){z). (A.48) 

n=oo 

It is therefore justified to write Qo{z) in the form / — - — . 

JT d,t 

3.4. Green's functions. 

Fix an element to in T and denote by Zq the space of functions z in L^'^ 
such that z{to) = 0. This is a (real) Hilbert space, with scalar product 

{zi\z2) = J^dtzi{t)z2{t), (A.49) 

hence Qo{z) is equal to | z). /,From (A. 44), one deduces the inequality 

\z{t)f<\t-to\-{z\z), (A.50) 

and therefore there is an element 5t in the dual Zq of Zq such that z{t) = {dt, z) 
(that is: 8t is a Dirac "function" centered at t). According to the general 
theory, one introduces a linear continuous and invertible map Gq : Tj'q ^ Zq 
corresponding to the quadratic form Qo on Zq. The Green's function is defined 
as follows: 

GQ{t,u) = {8uGoSu). (A.51) 

That is, the function t i— > Go{t,u) belongs to Zq and is equal to GqSu- By 
definition, one gets 

z{u) = {Go5u, z) (A.52) 

for every z in Z. More explicitly, these conditions on Go{t, u) can be expressed 
as follows: 

Go (to, U) =0 

j^dt—Go{t,u)—z{t) = z{u). (A.54) 
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The unique solution to these equations is given by 

{inf {t — to,u — to) for t >to, u> to, 
inf {to — t,to — u) for t < to, u < to, (A. 55) 

otherwise . 

We consider some special cases: 

- space Za of functions z in L^'^ with z{ta) = 0; 

- space Zb of functions z in L^'^ with z{tb) = 0; 

- space Za 5 = Za n Zfc with boundary conditions z{ta) — z{tb) — 0. 

The corresponding Green's functions are as follows: 

Ga{t, u) = e{t -u){u- to) + e{u -t){t- to) (A.56) 
Gi,{t, u) = d{t - u) {tb -t) + e{u - t) {tb - u) (A.57) 

Ga,b{t, U) = 9{t -U){t- tb) {tb - ta)~^ {ta - u) 

- 9{u -t){t- ta) {ta - tb)~^ {tb - u) . (A.58) 

Remark. Equation (A.58) can be obtained by integrating 

Vu;{z) S{z{ta)) z{t) z{u) . 
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3.5. Vector-valued functions. 

As in paragraph II. 1.1, we consider vector functions z — {z^,---,z'^) 
with components in L^'^. We introduce a real symmetric matrix (/Iq,/^) with 
an inverse denoted by {hf^^). The basic quadratic form is given by 

Qo{z) = [ dtKpz^{t)zl^{t). (A.59) 

The auxiliary condition 2; (to) = defines the space Zq. Similarly, the spaces 
Zo, Z5 and Za,b are described by the respective boundary conditions: 

- Z{ta) = for Za-, 

- z{tb) = for Zb-, 

- z{ta) = Z{tb) = for Za,b- 

In each case, the linear form 5" taking z into z'^{t) belongs to the dual of the 
corresponding of paths, and the Green's function is characterized by 

Gf{t,u) = {d^,Go6^) . (A.60) 
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Explicitly, we obtain 

Gf{t,u) = h'''^Go{t,u) (A.61) 

where the Green's function Gq refers to the boundary condition of the relevant 
space of paths (see formulas (A. 55) to (A. 58)). The Green's function can also 
be obtained from formula (A. 40) by specialization: 

/ Vz ■ e-'^^"^^)/^ z'^it) zf^{u) = ^Gf {t, u) . (A.62) 
Jz 27r 

Here Vz stands for Vg^Q^z. 



Remark. In general, the Green's functions G satisfying a second-order differ- 
ential equation DG = 1 are uniquely determined by d conditions at t = ta 
and d conditions ait = tj,. These conditions are obvious from equation (A.62) 
for the space Za,6, namely G°'^{t,u) = if t (or u) is one of the end points 

^a; ^fa- 
in the case of the space Z^, the boundary conditions are 

According to equation (A.62), these conditions can be expressed as 

/ Voj{z)-Z°'{ta)z^iu) = 0, 

I Vu;{z)- z''{tb)z^{u) = 0. 

The interpretation is as follows: any path z in satisfies z°'{ta) = 0; the 
time derivative i"(t6) at tb is totally unspecified, nevertheless it vanishes in a 
statistical sense being uncorrelated to the position and velocity at any other 
time. 

On the space Zq dual to Zq, we have defined the quadratic form Wq 
inverse to Qq. The Greens's function can be expressed as follows: 

Gf{t, «) = 1 [Wo{S^ + 6^) - WoiSr) - Wo{S^)] . (A.63) 

Conversely, given an element z' of Z' represented by 

{z',z) = I dtz'^{t)z'^{t), (A.64) 
Jtj 
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the quadratic form is given by 

Wo{z')= [ dt [ dtduGQ^{t,u)z'^{t)z'f3{u). (A.65) 
Jtj Jtj 

3.6. Scaling the paths. 

Here is a concise dimensional analysis of our quantities. 
Conventions : C (T) stands for the dimension of length (time) and [X] for 
the dimensional content of a quantity X. 

Since Qo{z) appears in e~^^°^'^^^^ , it has to be a pure number, hence we 
can use formula (A. 59) to deduce [hap] = C~^T. From formulas (A. 55) to 
(A. 58) we infer the dimension of Go{t, u) to be T. ^From (A. 62) we infer that 

Gq'^ = hence [/i"^] — C?T~^ from (A. 61). This is in accordance with 

the matrix relation haph^'^ = 5^ . Notice also that for a particle of mass m 
in a flat space, we have 

hai3 = mSap/h , h"'^ = hSa/3/m (A.66) 

(where h is the Planck constant and S^is the Kronecker delta), and that the 
dimension of h/m is C^T~^. We summarize our findings : 

TABLE 1 



Quantity 




t 




Go{t,u) 


Gf{t,u) h^f^ 


Dimension 


c 


T 


£-2r 


T 


£2 £2^-1 



We can confirm these results by scaling our paths. Let A > be a 
numerical constant, and denote by Az the path with components Az"{t) at 
time t. The basic integrator being written as T)uj{z) = e~^^°^^^/^'Dz, we scale 
it into an integrator 'Dio^{z) according to the formula 



Iz, 



Vu^{z)F{z) 



Vuj{z) F{Kz) 



(A.67) 



By suitably specializing F(2;), we obtain the new Green's function: 

/ Vu\z) z'^{t)z^{u) = A^Gfit, u) , (A.68) 
Jz 



in accordance with 



g: 



a/3 



C?. Similarly any n-point function scales as CJ 



The map taking z into Az is a linear map from Zq to Zq. By the general 
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theory in paragraph A. 2. 2, Vuj^ is another Gaussian integrator. From the 
formulas (A. 65) and (A. 68), it corresponds to the quadratic form A?Wq{z') 
on Zq, with an inverse quadratic form on Zq given by 

A-2go(^) = j dtA-^ ■ha/3z''{t)zf^{t). (A.69) 

This relation confirms [/ia/?] = A similar analysis applies to the scaling 

of time. 

We can rewrite formula (A. 67) in the following form: 

/ V^z-e-''^~'^°^'^/'H{A-h) = f Vz-e-^'^^^'^/'Hiz), (A.70) 

with a new integrator T>^z which is invariant under translations. It is justified 
to summarize the previous formula by V^z = V (A~^z), hence T>u)^{z) = 
Vlo(^A~^z) according to formula (A. 67). In the standard heuristic derivations, 
one writes Vz in the form C ■ Yit « dz^{t). According to our normalization 

I Vz ■ e-'^Qo(^)/* = 1, (A.71) 



we should write 



\[dz'^{t) 

t,a 

I J]d^"(t)-exp(-7rQo(-2)/s)' 



Vz = '-^ , (A.72) 



and similarly 



\{dz'^{t) 



t.a 



V^z = — . (A.73) 

/ W dz^{t) ■ exp (-ttQo (A"'-2) /s) 

Replacing z by A~^z amounts to replacing dz°'{t) by A~^dz°'{t), hence the 
volume element Ylta'^^'^i^) multiplied by A~^ where is the (infinite) 
number of degrees of freedom t, a. If we calculate V (A~^z) in accordance 
with (A.72), both numerator and denominator acquire a factor A~^ which 
drops out, and the correct formula V^z — V (A"^^;) is obtained from the 
heuristic formulas (A.72) and (A.73). 
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The heuristic constant A~ is equal to oo, 1 or according to the three 
cases 0<A<1,A = 1,A>1. This is reflected in the rigorous theory by the 
following fact : for K ^ 1 no functional F[z), except the constant 0, is such 
that both F{z) and F{Az) can be simultaneously integrated w.r.t. Vz In a 
finite-dimensional space IR^, the volume element dx^ ■ ■ ■ dx^ is the only one, 
up to a multiplicative constant, which is invariant under translations. In our 
infinite-dimensional setup there exist many translation-invariant integrators, 
but they act in different functional sectors. 

3.7. White noise representations. 

We consider now the Hilbert space L^(T) (or L^) consisting of the square- 
integrable functions ^ : T — > IR, with the standard quadratic form 



As usual we identify L^(T) with its own dual by means of the following scalar 
product 



Let us denote for a while the space L^(T) by H. Under our general conven- 
tions, both D : H ^ H' and G : H' ^ H are identity operators, hence the 
quadratic form H is its own inverse. The corresponding covariance is given 
by the kernel of the identity operator, namely 5{t — u). Hence denoting by 
the basic integrator, we get from equation (A. 40): 




(A.74) 




(A.75) 




(A.76) 



Let Lq'^(T) (or Lq^) be the subspace of functions z in L^'^(T) such that 
z{tQ) = 0. The spaces L^(T) and Lq'^{T) are isomorphic under a correspon- 
dance ^ 2; where ^ = z and conversely z{t) = J^^dt'^{t'). This can be 
expressed by the formulas 




(A.77) 




for 



for 



to<t' < t, 
t<t' < to, 



(A.78) 



otherwise. 



76 



The following transformation rule holds for functional integrals: 

/ Vz-e-''^°^'^/'F{z) = f V^-e-^'^^^^/'^i^), (A.79) 

where $(^) is equal to F{z) if ^ corresponds to z by (A. 77). The covariance 
in Lq'^ is obtained by specializing F{z) to z{t)z{u) in (A.79). Using (A. 76), 
we obtain 

Go{t,u) = J dt'e{t',t)e{t',u) (A.80) 

in agreement with formula (A. 55). 

To make contact with the heuristic definitions, we introduce the "coor- 
dinates" Xt = $,{t)y/di for the function ^, hence 



H{0 ^Y.^t , n = U (dXt/V^) . (A.81) 
t t 

We urge the reader to substantiate these claims by resorting to a subdivision 
T of the time interval T, like in paragraph A. 3. 3. 

Let X be any closed vector subspace of i^^(T), and let 11 be the orthogonal 
projector from L^(T) onto X, represented as an integral operator with kernel 
n(t, u). The quadratic form ffx on X is obtained by restriction of H, namely: 

H^{^) = J^dtm^ (A.82) 

for ^ in X. Denote by X'x'C the corresponding integrator on X. Then the 
covariance is expressed as follows: 



For instance, if we denote by -C^oC^) the subspace of L^(T) defined by the 
condition Jrj, dt^{t) = 0, the orthogonal projector is given by 

miu) = e(«) - -\- [ dt^t), (A.84) 

hence 

U{t, u) = S{t -u) — . (A.85) 

tb — ta 
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The map ^ ^ z where z{t) = dt' ^{t') takes isomorphically the space i^o(T) 
onto the space denoted by Za,b in paragraph A. 3. 4. It follows that the Green's 
function Ga,b corresponding to the space Z^^b is given by 

Ga,b{t,u)= f dv I dv'Il{v,v'). (A.86) 

Using (A. 85), we obtain easily (A. 58). 

Explicit formulas for transformations, mapping a given quadratic form on 
L^(T) into any given quadratic form on L^'^(T), can be found in [5, p. 274]. 
The reader is urged to extend the results of this paragraph to vector-valued 
functions. 



3.8. Fixing the endpoints. 

We denote by VaZ, V^z and Va^h^ respectively the integrators on the 
spaces Za, Z5 and Zo,6 (see paragraph A. 3. 5.). These integrators are related 
by the following formula: 

I Vaz ■ e-''^°^'^/'F{z)S{z{tb)) = C-^ I Va,b z ■ e-'^^°^^^/"F(z) , (A.87) 

where the constant C is given by 

C={detsG:'^{tb,tb)f' ^^^^^ 

A similar formula (where the roles of ta and tb are exchanged) can be found 
in [5, pp. 284, 357]. 



We give a new derivation of this formula to illustrate our methods. We 
simplify the notations by putting Vujaiz) equal to VaZ ■ e~'^^°^^''^^ , and defin- 
ing Vu)a,b{z) similarly. According to our general strategy, we need only to 
prove formula (A.87) in the special case 

F{z) = e-2"^("''") (A.89) 

where z' is an element of Z^. That is, we have to establish the following 
equality: 

/ Vu;a{z)-e-''''<'''')S{z{tb)) = C-^e-'''^'^''^''\ (A.90) 
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where Wa,b is the quadratic form on inverse of the restriction of Qq to 
the subspace Zo^;, of Z^. 

To prove (A. 90), we use the weU-known formula 

6{x) - / du . (A.91) 



We shall consider in detail the scalar case d — 1 and leave the general case to 
the reader. The left-hand side L of (A. 90) can be rewritten as 

L= du Vuaiz) exp (-27rz {z' + uSt„ z)) , (A.92) 

hence 

L= I du exp {-TTsWa {z' + u5tj) (A.93) 
m 

by (A. 18). We develop now the exponent 

Wa (z' + U5t,) = u'^Gaih, h) + '^uj^dt z'{t)Ga{h, t) + Wa{z') (A.94) 

and use the integration formula 

f du e-'^K^^'+s&^+c) = 1 exp i-TTs (c - h'^/a)) (A.95) 

to obtain 

L = ^ {AM) 
^ sGaitb, h) 

exp(^-TVs(Wa{z')-Ga{tb,hr^ J dtj du z' {t)z\u)Ga{h,t) Ga{tb,u)j^ . 

This can be transformed in the desired right-hand side of (A. 90) with C = 
■sJsGaith, tb) provided we establish the following identity: 

Ga,bit, U) = Gait, U) - Gait, tb) Gaitb, tb)'' Gaitb, u) . (A.97) 

We leave it as an exercice (see formulas (A. 56) and (A. 58)). 

We can prove a more general formula by a similar reasoning. We fix 
points Za and Zb and consider L^'-*^ paths z(-) with zita) = Za- We obtain the 
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afRne space + Z^, and transport to it the integrator T>aZ. Then, with a 
suitable constant Ca,b one gets 

I VaZ ■ e-"^°(^)/^F(^)5 {z{tb) - Zb) 

= C-\ \ T^a,bC ■ e-Q°(^)/^F (^ei + C) • (A.98) 

Here z^i is the afEne-hnear map Zc\{t) = Xt + fi with endpoints z{ta) = Za, 
z{tf)) = Zb, and the path ( in Za,bj that is C(^o) = Ci'^b) = 0, is a quantum 
fluctuation. 

4. Infinite-dimensional determinants. 

For this section, we refer the reader to [24], [25, chapter 5] and [26]. 

4.1. The case of operators. 

We return to our Banach space X. We assume that there exists on 
X an invertible^^ positive-definite quadratic form H (for instance H{z) = 

/t'^^^"(0^ the space Zq). We use the norm defined on X by ||a;|| = 
H{xy/^, and derive from it the dual norm ||a;'|| on X' as usual. 

According to Grothendieck, an operator T in X is called nuclear if it 
admits a representation in the form 



Tx^J2{x'^,x)-Xn (A.99) 

n>0 

with elements Xn in X and x'^ in X' such that J2n>o II^t^II ' ll^nll be finite. The 
g.l.b. of all such sums J2n \\-^n\\ " ll^^nll called the nuclear norm of T, to be 
denoted by \\T\\^. The nuclear operators in X form a Banach space, denoted 
by jC^(X.), with norm On C^CX), there exists a continuous linear form, 
the trace, such that 

TT{T) = J2K^^n) (A.lOO) 

n>0 



^ A positive-definite continuous quadratic form is not necessarily invert- 
ible. For instance, let X be the space £^ of sequences {xi,X2, • • •) of real num- 
bers with J2'^=i (^n)^ < oo and define the norm by ||a;||^ = (xn)'^- We 
can identify X with its dual X', the scalar product being given by ^^^^ x'^Xn. 
The quadratic form Q{x) = X^^i {xn/nf' corresponds to the map D : X — 
X' taking (xi, X2, • • •) into (xi/1, a;2/2, • • •). The inverse of D does not exist 
as a map from into £^ since the sequence 1, 2, 3, ... is unbounded. 
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for an operator T given by (A. 99). 



We introduce a power series in A, namely: 

'yp{T)y := exp (a IV(T) - ^ T,{T^) + ^ T,{T^) (A.lOl) 

By using Hadamard's inequality about determinants, we obtain the basic es- 
timate: 

\ap{T)\<f>/^\\T\\l/p\. (A.102) 

It follows that the power series 'Ylip>o'^vO^)^^ infinite radius of con- 

vergence. We can therefore define the determinant as follows: 

Det(l + T) := ^ c7p(T) (A. 103) 

p>0 

for any nuclear operator T. By homogeneity, we obtain more generally 

Det(l + AT) = ^ ap(T)AP . (A.104) 

p>0 

The fundamental property of determinants is, as expected, the multiplicative 
rule: 

Det (C/i o U2) = Det (C/i) Det {U2) , (A.105) 

where Ui is of the form 1 + T^, with Tj nuclear (for i = 1,2). From this, and 
the relation (Ji{T) — Tr(T), we get a variation formula (for U — 1 and SU 
nuclear): 



Det{U + dU) 
Det{U) 



1 + Tr {U-^ ■6U)+0 (\\SU\\fj . (A.106) 



Otherwise stated, if U{i') is an operator of the form 1 + T{i'), where T{i') 
is nuclear, depending smoothly on the parameter i/, we get the derivation 
formula: 

^inDet{U{iy)) = Tr(u{iy)-''-^U{iy)] . (A.107) 



Remark. For any other norm ||-||^ defining the topology of X, we have an 
estimate 

\\x\\ < \\xf <C\\x\\ , (A.108) 
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with a finite numerical constant C > 0. It follows easily that the previous 
definitions are independent of the choice of the particular norm ||a;|| = HixY^"^ 
in X. 

4.2. Explicit formulas. 

Introduce a basis {en)n>\ of X orthonormal for the quadratic form i?, 
hence Hi^^tnen) = Z^n^n- operator T in X has a matrix {tmn) such 
that 

Ten = J]e^-t^n. (A. 109) 

m 

Assume that T is nuclear. Then the series Ylm ^nn of diagonal terms in 
the matrix converges absolutely and the trace Tr(T) is equal to as it 

should be. Furthermore (Jp{T) is the sum of the series made of the principal 
minors of order p: 

(jp{T) = V det {ti^,i^) i<„<p . (A.llO) 

ii<---<ip 

For the operator U = 1 + T, with the matrix with elements Umn = Smn + tmn, 
we obtain the determinant as a limit of finite-size determinants: 

Det(C/) = lim det (umn) i<m<N . (A. Ill) 

N=00 l<n<N 



As a special case, suppose that the basic vectors are eigenvectors for 
T, namely 

Ten = Ken ■ (A.112) 

Then we get 

Tr{T) = J2K , Det(l + T) = J](l + An) 

n n 

where both the series and the infinite product converge absolutely. 

The nuclear norm \\T\\^ can also be computed as follows: there exists 
an orthonormal basis (cn) such that the vectors Tsn are mutually orthogonal 
(for the quadratic form H) and then ||T||j^ = llTenll- 

Remark. Let T be a continuous linear operator in X. Assume that the series 
of diagonal terms tnn converges absolutely for every orthonormal basis. 
Then T is nuclear. When T is symmetric and positive, it is enough to assume 
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that this statement holds for one given orthonormal basis, and then it holds 
for all. There are counterexamples when T is not symmetric and positive. 

4.3. The case of quadratic forms. 

Contrary to a widespread misbelief, there is no such thing like the de- 
terminant of a quadratic form. Consider for instance a quadratic form Q on 
some finite-dimensional space with coordinates x^, - ■ ■ ,x^, namely 

Q{x) = hc^px'^x^ . (A.113) 

If we introduce a new system of coordinates x-*^, • • • such that = w^x^, 
then we obtain 

Q{x) = hx^,x^xt' (A. 114) 

with a new matrix 

= W^u^^h^P ■ (A.115) 

The determinants D = det (/ia/j) and D = det {hx^j) are connected by the 
scaling relation: 

D = D-{dei{ul)f . (A.116) 
Hence, what makes sense is the ratio of determinants 

det(/^i^j)/det(/^gJ) (A.117) 
associated to two quadratic forms 

Qo{x)^hfpx'^xf' , Qi{x)^h^^lx'^x^ (A.118) 
on the same space. We shall denote it by det((5i/(5o)- 

To obtain an intrinsic definition, let us consider two continuous quadratic 
forms Qq and Qi on a Banach space X and assume that Qq is invertible. We 
associate to Qq and Q\ two continuous linear maps Dj : X ^ X' such that 

{DiXi,X2) = {DiX2,xi) , Qi{x) = {DiX,x) , (A.119) 

for X, x\ and X2 in X. By assumption, Dq is invertible, hence there exists a 
unique continuous operator [/ in X such that Di = Dq o U. We denote U by^° 
Qi/Qo- In case the determinant of U is defined, namely when L'^ — 1 is nuclear, 
we denote it by Det {Qi/Qo). Of course, when X is finite-dimensional, this 
definition agrees with the previous one and 

Det [Qi/Qo) = det (Qi/Qo) (A.120) 



^° A better notation should perhaps be Qo\Qi. 
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in this case. 

A procedure to calculate this determinant is as follows. Let V be a finite- 
dimensional subspace of X. By restricting Qq and Qi to V we obtain two 
quadratic forms Qoy and Qiy on V. Assume now that V runs through an 
increasing sequence of subspaces, whose union is dense in X, and that Q^y 
be invertible for every V . Then 

Det (Qi/Qo) =lini det {Qiy/Qoy) . (A.121) 

4.4. A pencil of quadratic forms. 
Here are our assumptions: 

- there exists an invertible positive- definite quadratic form on X; 

- Qq and Qi are continuous quadratic forms on X; 

- Qo is invertible; 

- if we set Q = Qi — Qo, the operator T = Q/Qq on yi is nuclear. 

The last condition is an intrinsic property of the quadratic form Q, 
namely the existence of a representation like 

Q{x) = J2^n{x'n,xf, (A.122) 

n>l 

where the series a„ converges absolutely, and the sequence of numbers 
((x^, x))^^-^ is bounded for every vector x in X. We express this property by 
saying that the quadratic form Q is nuclear. 

We say that A is an eigenvalue of a quadratic form Q' w.r.t. Qo if it is 
an eigenvalue of the operator Q'/Qq. This is tantamount to saying that the 
quadratic form Q' — XQq is not invertible. 

We interpolate between the quadratic forms Qo and Qi = Qo + Q hy 
putting Qi, = Qo + lyQ where v is a, real or complex parameter. The determi- 
nant 

A{u) = Det(Q^/Qo) (A.123) 

is defined, being equal to Det(l + uT). This is an entire function of the 
complex variable u, hence it vanishes for a discrete set of values of u (possibly 
empty). Furthermore, A is an eigenvalue of Q w.r.t. Qo if and only if the 
quadratic form Q-i/\ is non invertible, that is if and only if A(— 1/A) = 0. 

As a consequence, the function of real variable u A^u) has, at most, 
a finite number of zeroes z^^i), • ■ ■ , z/^p) in the interval [0, 1] with < z/(i) < 
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• ■ ■ < z/(p) < 1. According to the relation (A. 107), between two such zeroes, 
the foUowing differential equation holds: 

^enBet{Q^/Qo) = Ty{Q/Q,) (A.124) 

where = Qq + i^Q- 

4.5. Some integration formulas. 

With the previous notations, consider any f in [0, 1] distinct from I'ii)-, • • • , 
. Hence the quadratic form Qi, on X is invertible, with an inverse quadratic 
form on X'. There exists an integrator x on X characterized by the 
formula 

/ V^x ■ exp (^--Q„{x) - 27Ti {x', x)^ = exp (-tts W„{x')) . (A.125) 

Our claim is that x is proportional to Vqx, namely that there exists 
a constant such that x — I{y)VQX. More explicitly, we assert the 
formula 

I V^x ■ F{x) = I Vox ■ F{x) , (A.126) 
ix ix 

and that a functional F{-) on X is Feynman-integrable for V^x if and only 
if it is Feynman-integrable for Vqx. The constant I{u) can be obtained by 
putting F{x) — exp (— f (5i/(x)) into equation (A.126), hence 

I{i^)-^ = J Vox ■ exp {Qo{x) + iyQ{x)j^ . (A.127) 
First case s = 1: 

According to our conventions, the quadratic forms Qo and Qi are positive- 
definite and invertible. Since the quadratic form Q = Qi — Qq is nuclear, it 
follows from the spectral theory that Qo and Q can be simultaneously diago- 
nalized. Hence there exists a basis {en)n>i such that 

(oo \ oo 

Y^tnen) = 5^(tn)' (A.128) 
n=l / n=l 
/ oo \ oo 



\n—l / n=l 

^00 



with real constants An such that ^^=1 I'^^l ^ Since Qi = Qo + Q is 
positive-definite, we have 1 -|- > 0, hence 

(oo \ oo 
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is again positive-definite and invertible for v in [0, 1]. Hence the determinant 
of QvIQq is defined and 

Det(Q^/go) = n (1 + ^^^) > 0. (A.131) 

The main result is given by the following formula: 

/(z/) =Det(g,/go)^ • (A.132) 

The proof is obtained without difficulty using equations (A. 127) to (A.131) 
and the approximation formula: 

Vox-F{x) = \im / d^tF{tiei + htAre„). (A.133) 

Second case s = i: 

Here the basic formula is given by 

= |Det (Q./Qo)r^' ^ind(Q./Qo) (A.134) 

where the indexlnd {Qv/Qo) is the number of negative eigenvalues of Q^, w.r.t. 
Qq. To simplify the statements, we shall assume that Qq is positive-definite. 
When V runs over the interval [0, 1], this index remains constant except when 
V is crossing a singular value I'^k)-, where it experiences a jump. Hence the 
index Ind (Qu/Qo) is a sum of local contributions from the exceptional values 
^(1) ) ■ ■ ■ 7 ^(p) 7 ci phenomenon reminiscent of caustics. 

The easiest proof of formula (A.134) is obtained by following a strategy 
initiated by Nelson and Sheeks in [6]. It works in both cases s = 1 and s = i. 
Starting from equation (A. 127), we obtain by derivation 



A. 
dv 



(A.135) 



that is 

d 



J^V.x ■ exp (^-^Q,(x)) {~Q{x)) , 
enlM = - [ V<jjJx)Q(x). (A.136) 



dv 

Expanding Q{x), according to formula (A. 122) and using formula (A. 39) for 
the covariance, we obtain 



n=l 
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which can be transformed easily into 



(A.138) 



According to equation (A. 124), we conclude 

= |^£nDet(Q^/Qo) • (A.139) 

It remains to study the shift in phase when ly goes through an exceptional 
value z^(fe). 

Remark. To simplify matters, assume again that Qo is positive-definite. Let L 
be an invertible operator in X, of the form 1 + T where T is nuclear. Putting 
Qi — Qo ° L, it can be shown that the quadratic form Q — Qi — Qq is 
nuclear and that Det (Qi/Qo) is equal to Det(L)^. The quadratic from Qi is 
positive-definite, hence Ind{Qi/Qo) = 0. ^From our main formulas (A. 132) 
and (A. 134), we derive 

Vs,Q,{Lx) = |Det(L)| • Vs,Q,x . (A.140) 

Coming back to the general situation, where = Qo + J^Q, the reader 
will find in [5, p. 277-279] examples of operators L{u) such that Qi, = Qo o 
L(i>). In such a case, we obtain VyX — |DetL(z^)| • Vqx. 

Let us summarize the main results obtained in this paragraph: 

Let Qo and Qi be continuous quadratic forms on the Banach space X. 
We assume that Qo and Qi are invertible and that Qi — Qo is nuclear. 

(A) Assuming that both Qo and Qi be positive-definite, we obtain 

[ Vq^x ■ e-'^^i^^) = Det(gi/Qo)"'/' • (A.141) 

(B) For the the oscillatory integral, we obtain 

I Vi^Q^x-e'''^'^''^ = \Det{Qi/Qo)\~^''^i~^'"^^''^°^ (A.142) 
ix 

where the index Ind{Qi/Qo) counts the number of negative eigenvalues of Qi 
w.r.t. Qo- 

Equations (A.141) and {A.142) justify the basic formula (A. 126) 

/ V^x ■ F{x) = / Vox ■ F{x) 
Jx Jx 

with given by (A. 132) or (A. 134). 
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JAppendix B 
Functional determinants of Jacobi operators 



Many functional determinants of Jacobi operators have been computed 
in works on semiclassical expansions. They are conveniently expressed in 
terms of finite-dimensional determinants of Jacobi matrices. We give here 
the backbone of the method, and applications to Jacobi operators along a 
classical path characterized by d initial conditions (position or momentum 
at ta) and d final conditions (position or momentum at t^). Details, proofs, 
applications and generalizations of the equations presented here are scattered 
in the literature and we give a few selected references at the end of this 
appendix. A general presentation on the properties of Jacobi operators will 
be found in the Ph.D dissertation of John La Chapelle [University of Texas 
(Austin) Ph. D. expected in 1995]. 

1. Jacobi fields and Jacobi matrices. 

A Jacobi field is a vector field along a classical path obtained by variation 
through classical paths. A d x d Jacobi matrix is built up from Jacobi fields; 
each column consists of the components of a Jacobi field. We give explicit 
constructions when the action functional is S{x) = Jrj, dt L{x{t) , x{t)) . 

Let {xciin)} be a 2(i-parameter family of classical paths (critical paths of 
the functional S) with values in a manifold M*^ 

xM-.T^M'^ , T=[0,T]. 

We introduce the notations 

Xci(t;/x) := (a;ci(/x))(t) (B.l) 
Xc\{t\ fi) := dxci{t; tj)/dt (B.2) 

x'^,^^{t■,^i)■.= ^x,,{t■,^i)/^^i^. (b.3) 

We choose ji = (//^, • • • , fi^^) to stand for 2d initial conditions that characterize 
the classical path Xci{fi), assumed to be unique for the time being; namely, 

S'{x,i{ii)) = (B.4) 

has a unique solution for a given set fj, of parameters. 
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By varying successively the 2d initial conditions, one obtains 2d Jacobi 
fields {dxci/d/J,"'} for cu e {1, • • • , 2d}. It is easy to show that 

S"{xM)-^-0- (B.5) 

We assume for the time being that the quadratic form S"{xc\{p)) ■ is not 
degenerate for variations ^ which respect the boundary conditions specified 
by /X. Therefore the 2d Jacobi fields are linearly independent. 

There are several convenient basis for the 2(i-dimensional space of Jacobi 
fields. Choose as initial condition 

IJ,= {Xa,Pa) (B.6a) 

where 

Xa := xita) (B.6b) 

Pa:=dL/dx{t)\,^,^ . (B.6c) 
The corresponding Jacobi fields are 

j'f'it) = dxl, {t; Xa,Pa) /dpa,(3 (B.7) 

r^(t) . (B.8) 

Having introduced Pa as one of the initial conditions, the reader suspects that 
we shall also use"^^ 

k,^{t) = dpcl,, {t; Xa,Pa) /dpa,i3 (B.9) 
L/3{t)^dp,i^.{t;Xa,Pa)/dx^^. (B.IO) 

Let the Jacobi matrices J, K, K, L be made of the Jacobi fields j, k, k, £ re- 
spectively 

J"^(t,ta)=r^(t) (B.ll) 

K^p{t,ta) = k''l3it) (B.12) 
Kj{t,ta)=kJ{t) (B.13) 

La(3{t,ta) ^iap{t). (B.14) 

The 2d x 2d matrix constructed from the four d x d blocks J, K, K,L is a 
solution of the equation defined by the Jacobi operator in phase space. The 



The classical momentum , {t; Xa,Pa) is obtained by evaluating dL/dx{t) 
at the point x{t) = x^ {t; Xa,Pa)- 
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properties of the Jacobi matrices are numerous. We note only their properties 
ai t = ta, which can be read off from their definitions : 



J{ta,ta)^0 , L{ta,ta)^0 (B.15) 
K{ta,ta) = l , K{ta,ta) = l. (B.16) 

The matrices K and K are indeed the transposed of each other 

Ka^{t,ta) = K('a{ta,t). (B.17) 

Let M, N, N, P be the matrix inverses of J, K, K, L respectively. Hence we 
have 

J^''{t,ta)M^^{ta,t)^6^ (B.18) 

and three similar equations expressed in condensed form as follows 

KN = ]1, KN = 1, LP^l. (B.19) 

The matrices M, A^, iV, P are the hessians of the corresponding action func- 
tions (Van Vleck matrices): 

and similarly and P are the hessians of S {xci{ti,) , Pci{ta)) and (S(pci(^6),Pci(^a)) 
respectively. 

2. Jacobi operators and their Greens's functions. 

The Green's functions of the Jacobi operators in phase space can be 
expressed in terms of the Jacobi matrices J, K, K, L; they include the Green's 
functions of the Jacobi operators in configuration space. For the sake of 
brevity, but at the cost of elegance, we consider here only the Green's functions 
of the Jacobi operators in configuration space. 

In the previous paragraph, Jacobi fields were obtained by variation 
through classical paths. In this paragraph, we consider a one-parameter family 
of paths x{X) satisfying d fixed boundary conditions at ta and d fixed boundary 
conditions at t^, referred in brief as "a" and "6". The corresponding path 
space is denoted by Va,bM, hence 

x{X) eVa,bM CVM . (B.21) 
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We assume that for A = 0, the path Xc\ = x{0) is the unique critical point of 
the action functional S restricted to Va,bM. We set 



x{X, t) := dx{X, t)/dt , x'(A, t) := dx{\ t)/dX , x'(0, t) = ^(t). 
The second variation of the action functional gives the Jacobi operator at a^ci^ 

= J^dt (Li«,i^r W + L2a,ii3r{t)) ^"{t) (B.22) 

+ J^dt [L,a,2^Cit)+L2a,2f3C{t)) i^{t) 

where ^2^,2/3 is d'^L/dx'^{t)dx^{t) evaluated at x^, etc... 

If we integrate (B.22) by parts, we obtain the Jacobi operator in its 
differential garb (second order differential operator on the space of vector 
fields ^ along Xd together with boundary terms at ta and tb). It is often 
simpler to work with the quadratic form S"{xci) • ^ written in (B.22). We 
call functional Jacobi operator the kernel corresponding to the quadratic form 
S"{xc\) ■ CC> namely: 

1 5'S"{xd) ■ g 

2 5^^{s)6m =:JaM^ci,.,t). (B.23) 
Its inverse will be called its Green's function G** {xci,t,u), namely 

dt Japixch s, t)Gf^^{xcu t, u) = 5Z5{s - u) . (B.24) 

Provided the quadratic form S"{xci) ■ is non degenerate, the functional 
Jacobi operator has a unique inverse. We say S"{xci) is non degenerate if for 
any ^ in Tx^{Pa,hM there exists r] in this vector space with 

S"{x,{)-ir^^Q. (B.25) 

This equation says that there are no (nonzero) Jacobi field in the tangent 
space T^^^Va,bM to the space of paths with boundary conditions a and b. 

We list below the Green's functions of Jacobi operators at classical paths 
with different boundary conditions. A more abstract formula could encode 
all cases. Explicit formulas may be more useful for applications. 
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i) The classical path is characterized by Pci{ta) = Pa, Xci{tb) = Xb 

G{t,s)^e{s-t)K{t,ta)N{ta,tb)J{tb,s) 

^ ^ (B.26) 

- Q{t - s)J{t, tb)N{tb, ta)K{ta, S) . 

ii) The classical path is characterized by Xci{ta) = Xa, Pc\{tb) = Pb 

G{t, S) = - t)J{t, ta)N{ta, tb)k{tb, S) 

-e{t-s)K{t,tb)N{tb,ta)J{ta,s). 

iii) The classical path is characterized by Xci(to) = Xa, Xci{tb) = Xb 

G{t, S) = d{s - t)J{t, ta)M{ta, tb)J{tb, S) 

-d{t-s)J{t,tb)M{tb,ta)J{ta,s). 

iv) The classical path is characterized by Pci{ta) = Pa, Pci{tb) = Pb 

G{t, S) = e{s - t)K{t, ta)P{ta, tb)K{tb, S) 

~ ~ (B./y) 

-e{t-s)K{t,tb)P{tb,ta)K{ta,s). 



3. Semiclassical expansions. 

We have given examples of WKB approximations in two cases: 

i) The classical path of reference is characterized by initial momentum, 
final position (paragraph III. 2). 

ii) The classical path is characterized by initial position, final position 
(paragraph IV. 1). 

Two other cases are often needed: 

iii) The classical path is characterized by initial position and final mo- 
mentum. The transposition from case i) is straightforward. 

iv) The classical path is characterized by initial momentum and final 
momentum. The critical points of the action functional are degenerate when 
the classical system is constrained by a conservation law. We refer the reader 
to the reference [15] for this case. 

The computation of semiclassical expansions along the lines of paragraphs 
III. 2 and IV. 1, but with general action functionals S", brings in a second 
variation 

'5"(xei)-ee = Qo(0 + Q(0 (B.30) 
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where 

QoiO = j^dtL2a,2pt{t)i\t) . (B.31) 
Provided the Legendre matrix 

{t),x,x{t)) = d'^L/dx'^.dxl, (B.32) 

is invertible, the Gaussian integrator defined by the quadratic form(B.30) can 
be handled by the same techniques as the Gaussian integrator defined by the 
quadratic form 

j^dtKpt{t)if{t). (B.33) 

The contribution of the second variation S"{xci) • the semiclassical 
expansion of ^'(tfc, x^) is 

j exp (-^(go(e) + QH)) = Det (Qo/(Qo + g))'/' . (B.34) 

Then one uses the Green's functions (B. 26-29) to identify the ratio of Jacobi 
matrices whose determinant is equal to Det(Qo/(Qo + Q))- The results are 
identical to the results obtained by discretizing the functional determinants, 
and reported in [16]. 

4. References for Appendix B. 

The most complete summary, to date, of Jacobi fields in phase space, 
including degenerate critical points of the action, can be found in Appendix 
A of 

- C. DeWitt-Morette "Feynman path integrals" Acta Physica Austriaca, 
Suppl. XXVI, 101-170 (1984). 

For the proofs of the results presented in the present paper, and for 
properties of Jacobi matrices used in the present paper, see [2, 5, 15, 16] and 

- C. DeWitt-Morette and T.-R. Zhang "Feynman-Kac formula in phase 
space with application to coherent state transitions" Phys. Rev. D 28, 2517- 
2525 (1983). 

More on degenerate critical points can be foinid in 

- C. DeWitt-Morette, B. Nelson and T.-R. Zhang "Caustic problems in 
quantum mechanics with applications to scattering theory" Phys. Rev. D 
28, 2526-2546 (1983). 

- C. DeWitt-Morette and B.L. Nelson "Glories and other degenerate 
points of the action" Phys. Rev. D 29, 1663-1668 (1984). 

A short introduction to the Hamiltonian techniques underlying the above 
results can be found in 

- P. Cartier "Some fundamental techniques in the theory of integrable 
systems" in Lectures on Integrable Systems (O. Babelon, P. Cartier and Y. 
Kosmann-Schwarzbach edit.). World Scientific, Singapore (1994). 
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JAppendix C 

A new class of ordinary differential equations 



The purpose of this Appendix is to extend to the L^'^ case the famil- 
iar theorems about the existence and uniqueness of solutions of differential 
equations, and to describe the parametrization of paths in a curved space by 
means of paths in a flat space. 

1. Solutions of differential equations: the classical case. 

We shall follow the usual strategy, as described in any standard text- 
book, for instance Bourbaki [27]. More precisely, consider a domain Q in the 
Euclidean space M"^ and a vector field in O associating to every epoch t in the 
time interval T = [ta^^fe] and to every point x in Q a, velocity vector v{t,x) 
in M"'. We assume that v{t, x) is a continuous function of t and x. The two 
basic remarks are as follows: 

a) any trajectory of the vector field v can be prolongated as long as it does 
not reach the boundary of Q; 

h)J by the m,ean value theorem,, if the absolute velocity \v\ is bounded by 
a constant V along a given trajectory leading from Xa at time ta to xi, at time 
tb, then the mean velocity ^^J' is bounded by V . 

^Prom these remarks follows the following existence theorem (Peano): 
Assume that fl is a closed ball centered at Xa of radius L, and that 

\v{t,x)\<L/T (C.l) 

holds uniformly for t in T and x in Q, where T = tb — ta is the length of T. 
Then there exists a solution x : T ^ Q of the differential equation x = v{t, x) 
with the initial condition x{ta) = Xa- 

For the proof, we construct first an approximate solution by the Euler 
method. Wc select epochs ti, ■ ■ ■ , tN-i such that ta < < ■ • ■ < ^TV-i < tb 
and set to = ta , tN — tb- Then we define inductively points xq, xi, ■ ■ ■ , xn in 
Q by xq = Xa and 

Xi = Xi-i + (ti - f (tj-i, Xi-i) (C.2) 

for 1 < i < A^; the estimate (C.l) guarantees that the points xq, xi, ■ ■ ■ ,xn are 
in Q. We then interpolate linearly in each subinterval tj], and generate a 
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function xr : T — > f2 depending on the subdivision T — {to < ti < • • • < ^at) 
of T. By the mean value theorem and the estimate (C.l) we obtain 



\xr{t) -xr{t')\ < V\t-t'\ 



(C.3) 



for t,t' in T, where V — L/T. We then use AscoU's theorem: it asserts the 
existence of a sequence of subdivisions T{n) of T, whose mesh /S.{n) tends to 
0, such that xr{n)it) tends to a hmit x{t) uniformly for t in T. Then, one 
shows that xr(n) satisfies an approximate integral equation 



XT{n) {t) -Xa- ds v{s, x{s)) 



(C.4) 



where lim Sn — 0. By uniform convergence, the limit function x{-) satisfies 

n=oo 

the integral equation 



x{t) = Xa+ / ds v{s,x{s)) , 

Jtr. 



(C.5) 



fully equivalent to the differential equation 



x{t) = v{t, x{t)) 



(C.6) 



Both the uniqueness of a solution to the previous equation and the 
continuous dependence of the unique solution on the initial position 
tablished by an analysis of the stability of trajectories. Suppose given another 
trajectory y : T — > f2 with initial position = y{ta)', such a trajectory exists 
if ya is sufficiently close to Xa- Hence we have the integral relation 

y{t) = ya+ I dsv{s,y{s)). (C.7) 

Furthermore, assume that the velocity field satisfies a Lipschitz condition: 

\v{t,x) — v{t,x')\ < k\x — x'\ (C.8) 

for x' in and t in T, with a fixed constant > 0. Denote by 5{t) the 
distance between x{t) and y{t). From (C.5), (C.7) and (C.8) one derives 

\5{t)-5{ta)\ <k [ dsS{s). (C.9) 
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By using an iteration procedure reminiscent of Picard's method, one derives 
the foUowing estimate: 

S{t) < 5{ta)e^^'-''^'^ (C.IO) 

{^^Gronwall's lemma"). If = y^, then 6{ta) = 0, hence 5{t) = for all t and 
x{t) = y(t): uniqueness. Furthermore if 5{ta) — |j/a ~^a| tends to 0, then 
y{t) tends uniformly to x{t) for t in the finite interval T. 

The results obtained so far are local. To get a global existence theorem, 
assume that AT is a compact manifold, of dimension n, and that X is a time- 
dependent vector field on A?", namely X{t, x) belongs to T^N for t in T, and x 
in A^. Assume furthermore that x) is continuous in t, x, with continuous 
first derivatives in x. Choose a point Xa in N and a coordinate system around 
Xa- Then X is expressed in this coordinate system by a function f : T x Q — > 
R", where Q is a closed ball centered at Xa-, with radius L. Since a continuous 
function is bounded on any compact set, there exists a constant r > such 
that x)| < L/2t uniformly for t in T and x in Vt. /,Prom the local 
existence theorem, it follows that for any point Xq in the open ball centered 
at Xa with radius L/2, and any to in T, there exists a unique trajectory 
X : [to,to + T] — > N such that a; (to) = xq- By compacity of N, we can 
cover N by a. finite number of such balls. Defining a as the minimum of the 
life-times r, we conclude that the existence and uniqueness of a trajectory 
X : [to,to + a] —>■ N with x{to) = xq hold for every point xq in N. 

Consider now the time interval T = [ta, tb] and subdivide it into subinter- 
vals [ti, ti+i] of lengths < a (for z in {0, 1, ■ ■ ■ , A^}) where to = ta, tN+i = tb- 
So given any point Xa in A^, there exists a trajectory xo{t) (for to < t < ti) 
such that xo{to) — Xa- Set xi = xo{ti) and consider the trajectory xi{t) (for 
ti < t < t2) with initial point xi. By repeating this procedure, we construct 
step by step the required trajectory x :T ^ N such that x{ta) = Xa- 

Proceeding backwards rather than forwards, we conclude that given any 
point Xb in A^, there exists a unique solution x{t) = x{t; Xf,) of the differential 
equation x{t) = X{t,x{t)) with final position x{ti)) = x^. For fixed t, we 
define the transformation E(t) in A" taking xj, to x{t\Xb), hence 

x{t\Xb)=Xb-T.{t). (C.ll) 

2. Solutions of differential equations: the L^'^ case. 

We sketch the necessary modifications. A more detailed account will 
appear elsewhere. 

Consider an L^'^ path a; : T — > f2; its action is defined by 




(C.12) 
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<A, (C.13) 



The mean value theorem is replaced by the following estimate 

I _ 1^ 

tb — ta 

which follows easily from Cauchy-Schwarz inequality. 

We make the following assumptions about the velocity field v{t,x). 

(A) The function v{t,x) is jointly-measurable in {t,x), and continuous 
in X for a given t, and furthermore there exists a numerical function 
V -.T ^ [0,+oo[ such that 

\v{t,x)\<V{t) (C.14) 

holds for X inO. and t inT. 

(B) Denoting by A the action frj,dtV{t)'^, the radius L of the ball Q 
centered at Xa, and the length T of the time interval T obey the estimate 

A<L^/T. (C.15) 



Under these conditions, one proves an existence theorem for the differ- 
ential equation (C.6) with an L^'^ solution x : T O. As a first step, one 
replaces Euler's approximation (C.2) by the following approximate solution 

xr{t) ^ XT{ti-i) + I dsv{s,xr{ti-i)) (C.16) 
Jti 

for t in the subinterval [ti^i,ti]. The inequalities (C.13) to (C.15) guarantee 
that this trajectory xrit) remains in the closed ball Q. Furthermore, the 
velocity xrit) of this trajectory satisfies the estimate |a;r(^)| < V{t) for t in 
T, and by (C.13), we obtain the estimate 

\xrit)-xr{t')f <A\t-t'\ . (C.17) 

We can then invoke Ascoli's theorem, and find a sequence of approximate 
solutions x-r{n) converging uniformly on T towards an L^'^ function x. Each 
approximation satisfies an approximate integral equation, and in the 

limit the integral equation (C.5) is obtained using Lebesgue's dominated con- 
vergence. 

For the uniqueness, we need a Lipschitz condition of the type 

\v{t, x) - v{t, x')\ < k{t) \x - x'\ (C.18) 
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where the integral Jrj, dt k{t)'^ is finite. We use again a variant of Gronwall's 
lemma. 

The global results are obtained for a compact manifold N by the reason- 
ings used at the end of paragraph C.l. Only minor modifications are needed. 

3. Parametrization of paths. 

We consider a manifold N and d vector fields -^(i), ■ ■ ■ , -^(d) on N. We 
assume that they are of class and linearly independent at each point of 
A^. For X in N denote by the vector subspace of T^N generated by 
• ■ ■ , X((i)(x). The collection of these vector spaces is a subbundle"^^ H 
of the tangent bundle TN to A^. We consider also a real symmetric invertible 
matrix {ha(3) of size d x d, and define a field of quadratic forms hx on by 

hx{X^a){x),X^^){x))=haf3. (C.19) 

Fix a point Xb in A^. Wc denote by V^^N the set of L^'^ paths x -.T ^ N 
which satisfy the following conditions: 

(A) The endpoint x{tb) is equal to Xf,. 

(B) For each epoch t in T, the velocity vector x{t) lies in the subspace 
Hx{t) ofTx(t)N. 

We define the action of such a path by 

A{x) = J^dthx(t){x{t),x{t)) . (C.20) 

Since the vectors X(Q,)(x(t)) (for I < a < d) form a basis of Hx{t)-i and 
from the hypothesis that the path x is of class L^'^, we infer that there exist 
functions in L^(T) such that 

x{t) = X^^^{x{t))z^{t) . (C.21) 

Furthermore, the function is the derivative of a function z"^ in L^'^(T) 
normalized by z'^{tb) = 0. The vector function z = (2;^, • • ■ z^) is an element 
of the space denoted by Z5 in paragraph A. 3. 8. This construction associates 
to a path X in 'P^AT a path z in with conservation of action: 

A{x) = Qo{z), (C.22) 

where as usual Qo{z) is equal to Jrj,dthaf3Z°'{t)z^{t). 

Assume now that N is compact. By using the theory of L^'^ differential 
equations sketched in paragraph C.2, it can be shown that we can invert the 

The letter H stands for "horizontal" . 
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transformation x y-^ z. Given any z in Zf,, the differential equation (C.21) has 
a unique solution x in V^^^N, and we obtain a parametrization 

of a space of paths in a curved space N by a space of paths in a fiat space W^. 
For z in Z;,, we denote by x{t, z) the solution of the differential equation (C.21) 
with endpoint x{ti,,z) = xj,- We can also introduce a global transformation 
'E{t,z) : N ^ N taking Xf, into x{t,z). If necessary, we include tb in the 
notation and denote this transformation by S(T; z) or ta, z) in the case 
t = ta- The chain rule 

ta, ZbaJ) = ^{tb, tc] Zbc) ' S(tc, ta] Zca) (C.23) 

is a consequence of the uniqueness of the solution of the differential equation 
(C.21). Here Zba is z^ tc is an intermediate epoch, and the paths Zbc '■ [tc, h] 
and Zca '■ [ta, tc] are given by 

Zbc{t)^z{t) , Zca{t) ^ z{t) - z{tc) . (C.24) 

Remark. The more general differential equation 

x{t) = X^c.Mt))z^{t) + Y{x{t)) (C.25) 

can be handled in a similar way. In this case, we replace V^^N by the space 
V^^^N of L^'i paths x:T such that x{tb) = Xb and x{t) - Y{t) belong 
to Hj.(^t^ for every t in T. 

These constructions are related to the Cartan development map. Take for 
M a compact Riemannian manifold, and let be the corresponding bundle 
of orthonormal frames; it is a compact manifold. Fix Xb in M, and a frame 
Pb — {xb, Ub) at Xb. Then, by the Riemannian connection, there is defined a 
"horizontal" subspace T^^N of the tangent space Tp^^N, and the projection 
n : A?" — > M induces an identification of T^N with Tx^M. Since Ub defines 
an orthonormal basis of T^^M, we obtain a basis X(i)(pb), • • • , -^(d)(Pb) of 
T^^N . This construction is valid for every point in TV, hence we define vector 
fields • • • , on N. By using the previous construction, we get a 

parametrization 

P:Zb^Vf^N. 

But the paths in Vf^N are nothing else than the horizontal liftings of the L/^'^ 
paths in M. More precisely denote by Vx^M the set of L^'^ paths x : T — > M, 
such that x{tb) = Xb. Then, the projection 11 : M induces a mapping 

of P^iV into VxbM, and this map is a bijection. To conclude, we 
get a diagram 

Zb^V^^N^Vx^M 

and by composition a parametrization of Vxi,M by Z5. This is the Cartan 
development map for L"^'^ paths. The standard theory works for paths. 
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Figure captions 



Figure A.l: "Linear change of variables" 
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